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Abstract 

We propose some spectral type collocation methods well suited for the approxi¬ 
mation of fourth-order systems. Our model problem is the biharmonic equation, in 
one dimensional and in two dimensions when the boundary conditions axe periodic on 
one direction. It is proved that the standard Gauss-Lobatto nodes are not the best 
choice for the collocation points. Then, we propose a new set of nodes related to some 
generalized Gauss type quadrature formulas. We provide a complete anlaysis of these 
formulas including some new issues about the asymptotic behaviour of the weights and 
we apply these results to the analysis of the collocation method. 
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[ 

I 

L . Int rodu ction. 

Spectral methods are well suited for the approximation of the solution of elliptic or parabolic 
1 type equations. They are known to be very efficient for second-order problems and their 
| approximation properties are consistant with the infinite order observed through numerical 
1 experiments. We refer to the recent book [CHQZ] for a review of most of the current results. 
However, and quite surprisingly, they have not been much considered for the discretization of 
fourth-order systems. The regularity of the solution of such problems is not in question, since it 
is generally higher than for second-order ones. On the other hand, the spaces of discrete functions 

I 

j seem specially appropriate to approximate the solution of high order equations, since they consist 
I either of truncated trigonometric series or of high degree polynomials that are both indefinitely 
I differentiable. Hence, they are contained in Sobolev spaces of any order, which is not the case with 
I finite element spaces. Conforming discretizations can thus be worked out easily. In addition, it can 
; be noted that approximating a linear fourth-order equation with constant coefficients via a 
i Galerkin method using spaces of high degree polynomials with respect to each variable gives 
optimal results. Only a few collocation methods have been implemented up to now for this kind of 

j 

j problems [0][He], and no numerical analysis has been provided for them (we refer to the review 
paper [BM3] for a survey of the strategies). We do think that the corner stone of collocation 
I techniques is the choice of the collocation nodes, i.e., the finite set of points in which the equation 
| will be exactly satisfied. In spectral methods, as first suggested by D. GOTTLIEB [Go], these are 
i always built from the nodes of a Gauss type quadrature formula, for two reasons. First, the 
Lagrange interpolation operator associated with these nodes has very good approximation 
properties. Second, the quadrature formula allows for writing a variational formulation of the 
\ discrete problem; then, the Strang lemma provides an estimate to compare its solution with the 
I exact one. In parallel, if a problem is stated in a variational formulation and not in a strong form, 
using this quadrature formula leads to consistant discretizations. 

For second-order problems, two possibilities arise : choosing the collocation set from the 
| nodes of a Gauss formula or from a Gauss-Lobatto formula. The difference relies on the fact that 
! this set contains some points of the boundary of the domain in the second case and not in the first 
lone. Due to the Dirichlet boundary conditions, the second choice turns out to be more efficient and, 


I 
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in two dimensions, only Gauss-Lobatto points lead to an optimal approximation error. 

In order to discretize fourth-order equations, the question is not so easy to solve since two 
boundary conditions, one on the function and one on its normal derivative, must be enforced at each 
boundary node. That is why we propose to use a generalized Gauss type quadrature formula, which 
approximates the integral of a function on a real interval by a sum of its values at some interior 
nodes plus its values at the extremities of the interval plus the values of its derivative at the 
extremities, each of them being multiplied by an appropriate weight. Note that the nodes of the 
Gauss formula are the zeros of a fixed orthogonal polynomial, those of the Gauss-Lobatto formula 
are the extrema of this polynomial, i.e., the zeros of its first derivative; by similar arguments, it 
turns out that the nodes of the generalized formulas must be chosen as the zeros of the successive 
derivatives of this polynomial. We shall thoroughly study these quadrature formulas, both from 
theoretical and numerical points of view. 

I 

Our aim is of course to discretize fourth-order problems by collocation techniques involving 
the nodes of the generalized Gauss type formula. We first consider the simple case of a j 
fourth-order equation on a finite real interval, when the solution must vanish at the extremities 
of the interval, together with its derivative. On this test problem, we compare two discrete 
problems: in both of them, the exact solution is approximated by a polynomial of the same degree 
which satisfies the boundary conditions, but the equation is enforced at the interior nodes of either i 
a Gauss-Lobatto formula or an appropriate generalized Gauss type formula. Finally, we consider , 

I 

the equation of the bilaplacian on a rectangle, when the boundary conditions are periodic in one 
direction and homogeneous in the other one. We discretize this equation by a collocation method 

j 

using the nodes of the generalized formula in the nonperiodic direction, and we provide a complete 

I 

numerical analysis of this method. Our theoretical justifications are all given in the generalized^ 
framework of weighted Sobolev spaces [BM2], which allows for a simultaneous treatment of the ( 
Legendre and Chebyshev collocation techniques (we refer to [CHQZ] for a comparison between 
them). Our intention is to extend our method to fourth-order problems in a rectangle, provided 1 
with inhomogeneous Dirichlet conditions, in a forthcoming paper . j 

An outline of the paper is as follows. Section II is devoted to the analysis of the generalized' 
quadrature formula. A variational formulation of the monodimensional fourth-order problem foH 

> 
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the bilaplacian in weighted Sobolev spaces is studied in Section III. In Section IV, we compare two 
collocation techniques for approximating this problem. Finally, in Section V, we extend the method 
to the two-dimensional equation with mixed periodic-nonperiodic boundary conditions. The paper 
contains three appendices: the first one gives general approximation properties of high-degree 
polynomials in the weighted Sobolev spaces of order 2; the second one states the approximation 
properties of the Lagrange interpolation operator at the nodes of the generalized Gauss type 
formula; the third one contains several tables of nodes and weights of the quadrature formulas. 
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Jl. - Th e generalize d , q uadrature fo r m ula, 


In all that follows, we denote by A the open interval ]-!,![. For any integer n > 0, P n (A) is 
the space of the restrictions to A of all polynomials of degree ^ n. 

For any real number a > -1, we define the weight p K on A by 

(11.1) VC€A, qJX) = ( 1 -C 2 )* • 

With this weight, we associate the following scalar product, which is defined on the space of all 
functions, the square of which is integrable with respect to the measure dC , 

(11.2) (•P.4>) k = ^, *p(t) *U) P*(Odt; . 


We recall that a family of orthogonal polynomials with respect to the scalar product (...^ is the 
family of Jacobi polynomials ( J* ) ncW . where J* has degree n and satisfies the condition 

01.3) j»(ii) = uo" . 

n! r(oc+ 1) 

(r denotes the classical Euler's gamma-function). A number of properties of these polynomials are 
well-known (see [DR, §1.13] or [Ho]): among them, we have 


(11.4) V (m,n) € IN 2 , J_, J“(t) J* (O P a (U d£ = 6, 


2 2ct+1 r(n+oc+1) 2 


mn 


(2n + 2oc + 1) n! r(n + 2a + 1) 


where 8 mn stands for the Kronecker’s symbol; the family (J* ) nc(N satisfies the induction formula 
(11.5) (nf 1)(n4 2a • 1) J* +1 = (2n.2a +1 )(nta+ 1) K J* - (n .oc)(n-fOc +1) J*_, , 

Jj(C)=1 et J*(t) = (oc +1) t; 

We shall also need the following result (see [BM2, Lemma IV.2]), which is valid for n > 2 


(11.6) Jj n *(t)dt = 


[ 


n + 2oc+ 1 


nm 


'n + 1 




2n*2oc«1 * nfoc+l n+l n + 2a 
where j J“(0 dC denotes the primitive function of J* which is orthogonal to 1 for the scalar 
product (...) . 


Moreover, each polynomial J* , n € IN, is an eigenfunction of the operator A a defined by 

(11.7) A tt v = - P_«. (P tt+ i v’)' , 

indeed, it satisfies the following ordinary differential equation 

(11.8) (p K+1 v) n Kl ) > + n(n + 2« + l)p K J* = 0 . 

From this equation, we observe a property which is the corner stone of our analysis : the family of 
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polynomials ( J*' ) , is orthogonal with respect to the scalar product (.,.) K+1 . hence, for any 

n ^ 1, the polynomials J*' and J**] coincide up to a multiplicative constant. More precisely, 
using (11.8) and (11.3) to compute J*'( 1), we have for any n ^ 1 

(11.9) J*’ = . n+2« + L . 

Now, let N be a fixed integer. In all that follows, we denote by , 1 < j < N, the zeros of the 
polynomial J(j (we drop out the index N for sake of simplicity); it is well-known that these zeros 
ore distinct, so we may assume that : < ••• < • With each zero , 1 < j < N, we 

associate a characteristic polynomial Q* , i.e. the only polynomial in P N-1 (A) which is equal to 1 
in and vanishes in , 1 i < N, i * j. For any function $ defined on A and any pair (p_ , p + ) 
of real numbers, the notation <P(± 1) p ± in a summation denotes the sum $(-1) p_ +#( +1) p + . 

We are interested in quadrature formulas to approximate the integral j_, $(0 P K (C) dC, 
where<£ is a function defined on A and is assumed to be smooth; moreover, we want these formulas 
to be precise, i.e. to be exact on polynomials of the highest possible degree. Two formulas are 
well-known: 

1) the Gauss formula 

f! 1 9(0 Q a (K) dt - Lj H , i *<£*) e*' 0 ; 

the nodes are the zeros of Jjj ; for a suitable choice of positive weights p*' G , the formula is exact 
on P 2 N _i<A). 

2) the Gauss-Lobatto formula 

1’, -mo 9«(odt =. z:" = ,wc;'•’> *fuo . 

A 

(with the sommation convention described above); here, the nodes are the zeros of (1 — C ) J N ^, , 
i.e. by (11.9) the zeros of j£ +1 and the bounds ± 1 of the interval; for a suitable choice of positive 
weights p“ ,GL , 1 ^ j < N, and p“ ,GL , the formula is exact on P 2N+1 (A). 

We propose the following generalized quadrature formulas, depending on a nonnegative 
integer m : 

(11.10) (', MO 9.(0 dC - L, N s , Mt J - m ) tr * E k ".'o (<^/dC l )(±1) p“;™ . 

(with the same sommation convention). Here, the interior nodes are the zeros of the polynomial 
j“ +m . Clearly, with a suitable choice of the weights, the Gauss formula is obtained for m a 0 and 
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the Gauss-Lobatto formula is obtained for m = 1. 

Next, we define the weights in order to make the formula as accurate as possible. 

Lemma LLi : For any real number cx > -1 and for any integer m > 0, there exists a unique 
(N-t 2m)-uple of weights p“' m , 1 < j ^ N, and p“'™ , 0 < k ^ m-1, in (R N x (R 2 " 1 such that the 
quadrature formula (11.10) is exact on P N+ 2 m _i(A). 

Proof: The quadrature formula (11.10) is exact on P N+ 2 m-/if and only if if is exact on a basis 

of P N + 2 m-i (A) - i- e ' the vector ( g“' m ,.... p^'" 1 . P?.'* . pJIm .* > is a solution (X, .X N , y 0± , 

.... y m _, 1 ) of the linear system 

Lj N : , «r> n Xj + L k m :d (d k (L n )/dL k )(± 1) y k<4 = f0 P*(0 dL , 0 < n < N + 2m-1 . 

This is a system of N + 2m equations with N + 2m unknowns, hence it has a (unique) solution if and 
only if the only solution of the same system with a zero right-hand member is zero. Therefore, let 
(X, ,.... X N , y 0± ,y m _ 1± ) be a solution of 

(11.11) V * € P N+2m _i(A), Lj N s , Xj . L k m :d (d k< I , /dC k )(± 1) y k± = 0 . 

First, choosing $ = (1 -L 2 ) m Q* +m in (II. 1 1), we obtain X. = 0, 1 < 1 < N. Next, for t? decreasing 
from m-1 to 0, we choose successively <P = (1 -L 2 ) 1 (1 +0 ar, d$ = (I -C 2 )* (1 -O in 

(11.11) , so that we deduce y, _ = y, + =0. That ends the proof. 

In all that follows, we assume that the weights p* ,m , 1 < j < N, and p*'™ , 0 ^ k ^ m-1, are 
chosen such that the formula (11.10) is exact on P N+ 2 m _i(A). We derive the 
Proposition LLI : For any real number <x > -1 and for any integer m > 0, the quadrature 
formula (11.10) is exact on P 2 N + 2 m-/ A ^ 

Proof : Let $ be a polynomial in P 2 N + 2 m-i^ A ^ &y the Euclidean algorithm, we can find a 
polynomial Q in P N _,(A) and a polynomial R in P N+2m -i(A) such that 
<P = (1-C 2 ) m Jn +R1 q + R 
Next, we compute 

(d, <KU P*(0 dc = id, (1 -L 2 ) m JN tm (0 0(0 P«(0 dc + j], R(0 P a (0 dC 
Since J^ +m is orthogonal to any polynomial of degree < N-l with respect to the scalar product 
(■ >• ) K+m ■ we obtain 
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W) e.(OcR; =1^(0 0.(0 dC . 

so that, by Lemma 11.1, 

/’,*«:) e,(U dt = L," t , R(c; ,m ) e“ m * L k m ;o tfR/dJXi I) p k *;” 

= l -, , ■Hep) pj , ' m * tail (d^/dt‘)d i ) pj” 

Remark 11.1 : Of course, the property of Proposition 11.1 is satisfied by the Gauss formula 
(m = 0) and by the Gauss-Lobatto formula (m = 1). The idea of building quadrature formulas 
involving some values of the derivatives of the function is not new (see [DR, § 2.7][St][T]). Note 
that these formulas are essentially of Gauss type since the choice of the N free nodes on (R that 
generate the quadrature formula is optimal: any other choice leads to a method which is not exact 
for all polynomials of P 2N+2m .i^)' However formulas of type (11.10), the nodes of which are 
zeros of Jacobi polynomials, are specially interesting since they can be used to introduce and 
justify new spectral collocation methods, as it will appear later. 


Since the nodes of formula (II. 10) are the zeros of a Jacobi polynomial, they can be computed 
by the same algorithms as for the Gauss formula [DR, § 2.7] : for instance, the zeros t;" , 
1 ^ j ^ N, are the eigenvalues of a tridiagonal symmetric matrix of order N. Details about this 
method and tables giving the values of these nodes for different values of the parameters a, m and 
N can be found in Appendix C. 

We end this section by giving some properties of the weights p* ,m , 1 ^ j < N, and p*'™ , 
0 k m-1. Our purpose is double : theoretical (for instance, the positivity of some of these 
weights will be very useful in several proofs) and numerical (since we need them in a number of 
practical computations). 


We begin with the "internal" weights p*' , 1 < j < N. Let us recall ([DR, § 2.7] or [Sz, 

formula (3.4.7)]) that the Gauss weights are given by 


(n.12) pr-° = 


K,0 / V N " ' 


«n2 




j n “(U 2 e.< Odt 


) 




Due to the Christoffel-Darboux formula [Ho, 12.1 ][$z, Thm 3.2.2], this is equivalent to 



! 
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(11.13) 


« o 

<v =- 


2 2k+1 r(N mx + 1) r(N + oc + 2) 


(N+1)! r(N + 2a + 2) 


> 2 * +1 r(N + a) r(Nm+ 1 ) 


N! r(N + 2«4 1) 


j;.1(cpj n k up 


The following lemma shows how to compute the weights p* ,m from p* +m,u , 1 ^ j < N. 


a-i-m.O 


Lemma 11.2 : For any real number oc > -1 and for any integer m ^ 1 , the weights p* ,m , 

1 j < N, satisfy 

(11.14) p* ,m □ (1 -(L“ +rn ) 2 )' m p* +m, ° 

Proof : For any $ in P N _,(A), let us compute the integral J 1 , (1 -t; 2 ) m #(0 P K (0 dt; both by 
formula (11.10) and by the Gauss formula associated with the weight p K+m . We obtain 

J ’, (1 -C 2 ) m <KO P a (U dC = Z* s , (1 -(^ +m ) 2 ) m p*’ m = Lj H a , #U“ +m ) p“ +m ' 0 . 

Choosing # = Q“ +m , we prove (11.14). 

Remark 11.2 : Formula (11.14) shows that the weights p“ ,m , 1 < j < N, retain some properties 
of the Gauss weights : indeed, for 1 < j < N, we have 
(ms) p;- m > o , 
due) p n -j= 


Next, we consider the "boundary" weights p“’™ , 0 < k < m- 1. Let us first compare p*’™ and 

0k.- 

Lemma 11.3 : For any real number <x > -1 and for any integer m > 1 , the weights p“'™, 

0 k ^ m-1, satisfy 
(11.17) p*; + m = (-1) k p*;™ 

Proof : Let# be any polynomial in P 2 N + 2 m-i^^ Using (II. 10) to compute \\ s #(L) p tt (0 dC = 
#(-C) P a (0 dC, we have 

C“ ,m ) e*' m . L?:'o (d^/dc'K-O (C . L k m :» (d k d>/dc l X + 1 ) p k ;™ 

= L-,, p;- m * <- > ) k \ > p k *-: : 

* LjT.o (-') k (df'*>/dC k )(—1) P k ” i 


l 



-9 


whence, bg (11.16), 

z?:i (d^/dt k )(-D <>:? * xr:d (dVdex.n 

= L k m -0 (-')‘ CdVdt. k )< 41) p k “;. m * zrj (-' ) l <df*/dt k )<-1) p™ 

For P decreasing from m- I to 0, we choose <1X0 = (1-t) m (1 tf;)* , which gives (11.17). 

The following lemma precises the sign of p£'™ 

Lemma 11.4 : For any real number oc > — 1 and for any integer m > 1 , the weights p£.'™ 
satisfy s 

(11.18) p“;_ m □ p 0 % m > 0 

Proof: Applying the formula (II. 10) to the function*? = 1 gives 

r 1 / y \ ,y V N a,m A a,m 

J_ 1 - 2-.j : i Pj + Po,± 

so that, bg Lemmas 11.2 and 11.3, 

2 e ; : . m = 2 ej? = L,o -c 2 r m 9 ..JO dc-X" s , (I -(c;* m ) 2 )- m 

We have proved that p£'™ = p£'™ is equal to the half of the quadrature error of the Gauss formula 
with N points for the weight p tt+m , applied to the function (1 -C 2 )~ m Recall the general result 
giving the quadrature error for the Gauss formula [CM, Th. 2.5] : for ang function <? in 0 2N ( A), 
there exists t, -1 < L < 1, such that 

|',f<0 0,(0 dt - , <Htp p“-° = (d 2N <P/dt 2N )(f;) IM;IIo,,. a / <2N)! k N « 2 , 

where k“ denotes the coefficient of t N in Jjj • Hence, to obtain (11.18), it suffices to prove that the 
2N- th derivative of (1 -t; 2 )"" 1 is nonnegative on A. Noting that 

(i -er m = ( l; = °°o c 2 * r . 

we see that this function and its 2N-th derivative are even and that all the coefficients of their 
series expansions are nonnegative. Consequentlg, theg are ^ 0 on A. 


! It is alreadg known [BM2, Lemma V.3] that, for the Gauss-Lobatto formula, the weights pg' ± 

are given bg 

{ (11.19) q^_ = qI‘1= l]iJS +, (0(W0p K (t)dC/2J* +, (1) , 

! which gields bg (11.3) and (11,6) (see [BM2, Lemma V.3]) 


( 11 . 20 ) pj ;’ = p *;’ = 2 2cu1 r ( a + 1 ) r(<x + 2 ) 


T(N + 2 a + 3 ) 
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The next lemmas show how to compute the weights p“'™ , 0 < k < m-1. 

Lemma 11.5 : For any real number <x > -1 and for any integer m ^ 1 , the weights p*'™, 

1 ^ k < m-1, satisfy 

( 11 . 21 ) 2 (m- 1 )o m :{ n * = -e^ : r 1 

and 

(11.22) 2k p“‘™ + k(k+1) p k *‘ 1 m + = - , 1 ^ k ^ m-2 . 

Proof: Let us choose^ = (1 -L 2 ) ^, where ^ is an even polynomial in We com P u ^ e 

j], <p(0 P K (0 dt; = J., ^(O p a+1 (C) dt; by using (II. 10) once for a and m, once for a + 1 and 
m- 1. We obtain 

L"., (I -(t;*"’) 2 ) Pj*” . z::l I -t 2 ) >K)/dc l ](± 1) 0 k ”;" 

= Lj s , W,;* m ) . L k m : 2 1) p k -’”-' 

Using (II. M), (II. 17) and the fact that ^ is even, we derive 

z":\ i*<' -c 2 > w^k. i) p k ;, m = l::1 <dVdc k x > i> . 

But it turns out that, for k ^ 2, 

[d k (( 1 -C 2 ) ^)/dt; k ]( + 1) = - 2k (d k_1 ^/dt; k " 1 )(+1) - k(k-1) d k " 2 ^/dt; k " 2 ( +1) , 

and that ((1 -t; 2 ) ^)'(+ 1) is equal to - 2 SL'(+ 1). Hence, we have 

- 2 ¥(+1) p“7 - t 2k (d k_1 ^/dC; k " 1 )( +1) + k(k-1) d k " 2 ^/dC k ' 2 ( + 1)] p*;™ 

=L:: 0 2 (^/dc k )( + i)e;:: jn - 1 

For F decreasing from m-2 to 0, we choose ^(O = (1 - t; 2 )* , which yields (11.21) and (11.22). 

Remark 11.5 : It follows from formulas (11.21) and (11.22) that the vector of weights p*’™ , 

1 < k < m- 1, can be computed from the vector of weights p* + ^' rn_1 , 0 ^ k < m-2, by solving a 
linear system. The matrix of this system is upper triangular and bidiagonal. 

Remark 11.4 : Due to formulas (11.1 7), (11.18), (11.21) and (11.22), we observe that, for 
0 < k ^ m- 1, 

(11.23) p«;™ = (-1) k p k “T > 0 , 

hence the only negative weights are the p“'™ for odd values of k. 
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The computation of the weights p£'™ involves a quite technical result. Let us set, for any 
a > - 1 and any integers N ^ 0 and m 1, 

(11.24) (d m J N K +m /dC m )(Up„(OdC . 

Of course, this quantity is equal to 0 when N is odd. But we are interested first in even values of N. 

Lemma 11.6 : For any real number a > -1 and for any integer m ^ 1 , for any nonnegative even 

integer N , the quantity X£ ,m is given by 

m 2a . m r(a+l) r(N + a+m +1) (N/2 + m-1)! r(N/2+a + m +1/2) 

Ml ?F) X ' - o 2a+m __ 

N (m-1)! r(N +2a + m +1) (N/2)! r(N/2 + a + 3/2) 


Proof : From formula (11.6), we induce 


J 


a 

N +m 


N t a t m 

(2N t 2a + 2m-1)- 

N 4 2a4m 


J 


a. 

N +m- 1 


(N+a »m)(N +a *m- 1) 

__ vj ** * 

(N + 2anm)(N+2a + m-l) N_2+m 


First, derivating this equation (m-1) times, we have at once 

N+a + m (N+a + m)(N+a + m-1) 


X“' m □ (2N + 2a + 2m- 1) 


N+2a+m 


y«,m- 1 

X N f 


(N + 2a+m)(N + 2a + m-1) 


wOt, m 
A N-2 


On the other hand, using the same formula with m = 1, we obtain 


v «.i _ (N + a^1)(N + a) ...(a42) yK>1 

N = (N + 2a+ l)(N42a) ...(2a + 2) 0 ’ 

= (r(N+a + 2)/r(a + 2)) (r(2a + 2)/r(N + 2a + 2)) (a + 1) j], p tt (0 dt; 
which implies by (11.4) 

X"' 1 = 2 2k+1 r(aH) r(N«at2)/r(N.2ai2) . 

Finally, using (11.9), we note that 

X“' m = (d m J“/dL m )(t;) P tt (0 dL = (2a + m 1 1 )(2a + m 1 2) ... (2a + 2m) Q*(0 dL/2 m 

= 2 2 “' m+1 (T(a+1 ) 2 /T(2a f2)) (r(2a + 2m f 1 )/r(2a + m + 1)) . 

The formula r(2s) = (1 /V2n) 2 2, ' 1/2 r(s) r(s* 1 /2) then implies 

I 

x «,m n 2 2K+m (r(a+ 1)/r(a + 3/2)) (r(a + m+1 )r(a+m +1/2)/r(2a + m +1)) . 

The induction formula on XjJ ,m , together with the values of Xjj' 1 and x£ ,m , allows us to prove 
(11.25). Indeed, to make the computation readable, let us set 

YjJ/J* = (r(N + 2a-fm + 1)/r(N+a + m+1)) X*' m / 2 2k+1 T(a+1) . 

The sequence ( Y£ ,m ) , K?0 satisfies 
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Y“ m = (4K + 2oc + 2m-1) Y*' m_1 + , m>2,K>1, 

Y"' 1 = I , Ks*0, and Y“- m = 2 m ' 1 r(oom + 1 /2)/ r(oo3/2) , m»1 . 

Then, it is easy to check (but not so easy to find) the formula, valid for any m ^ 1 and K > 0, 

Y«,m _ 2 m-i (K^-0! r(Ki(X *m + 1 /2) 

K K! (m-1)! f(K »cx+3/2) 

which yields (11.25). 


An immediate consequence is the 

Corollaru 11.1 : For any real number a > -1 , for any nonnegative even integer N , the 
sequence ( X^ ,m ) m , is given by 


(11.26) 


X*' m - (N4 2onm * 1 )(N + 2m-2)(N+2co3) X 


cc+ 1 ,m- 1 
N 


/ 8 (oc + 1 )(m-1) , m 2, 


/«,i 


-,2a* 1 


r(a . 1) r(N+a + 2) / T(Nf2oo2) 


We also need the 

Lemma 11.7 : For any real number a > -1 and for any integer k ^ 1 , for any nonnegative even 
integer N , the polynomial Jjj satisfies 

(11.27) (d k J*/dt k )(+D = ((N +1 -k)(N + 2cx + k)/2(ook)) (d k ' , J l ^/dt; k ‘ 1 )( + 1) . 


Proof : Formula (11.8) gives 

(1 -C 2 ) vJ“" - 2(00 1) LJjJ' + N(N + 2ool) = 0 . 

Then, it can easily be proven by induction 

(1-C 2 ) (d k+1 J^/dt, k+, ) - 2(a»k) Ud k J*/dt; k ) + (N+1 -k)(N .2a *k) (d k “ 1 J*/dt; k_ 1 ) = 0 
whence the result. 


Formula (11.27) leads us to define, for any oc > -1 and any integer N > 0, the sequence 

0 by 


(11.28) 


X t * = (N+ 1-k)(Nf2ook) X k a _, / 2(cok) , k>l, 


x;= i 
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I fttnma 11.8 : For any real number a > - 1 and for any integer , the weights are 

given by 

1) if N is even, 

/II ->n\ „<X,fTl -«.m 7*.^ V m ”1 \ a + m n K ' m 

(11.29) P 0 ,_=P 0 . + = Z -L k= i\ p k , + 
where the sequence ( Z a,m ) m ? , is defined by 

: «,m z (N.2m-2)(Nt2cc*3) _ m > 2> 

(11.30) 4 (a + 1 )(m- 1) 

Z*' 1 = 2 2a+1 r(a+1) r(a + 2) N! /r(N+2a + 3) ; 

2) if N is odd, 

/ii „a,m T <x,m v ( \ * +fn » k \* + m 'l n <lt ' m 

(11.31) p 0 _ = P 0j+ = T -2* k = 1 ( A k +KA k -l ; Pk, + 

where the sequence ( T a,m ) m „ , is defined by 

(N*2m-3)(N.2o<.4) 

j / I 4 \ 1 


01.32) 


4 (a +1 )(m-1) 

T a, i = 2 2a+1 r(a+1) r(oc + 2) N!/r(N+2a+3) ; 


Proof : First note that, due to Lemma 11.7, we have 

X k K+m = (d m+k J^ m /dC m+k )( +1) / (d m J“ +m /dC m )( + 0 ■ 

1) When N is even, applying formula (11.10) to the polynomial <P = d m j£ +m /dt; m gives at once 

x;*' = 2 <cTj M “, m /dc m K. i) . 2 l ::1 (<r ,k j; im /dc m *")( . t) eC:r 

The result follows by setting Z*' m = x£' m /2(d m JN +m /dt; m )( +1) and using Corollary 11.1 and (11.9). 

2) When N is odd, we must apply formula (11.10) to the polynomial <? = L d m J|J +rr /dt; m We obtain 

j’, t (d m j N *. m /dr)(o p.(t) * = 2 (<rj;. m /dc m )(.t> Po * : : 

, 21";.' {(d m ‘ k o;. m /dt; m * k )(*l) ♦k(d m,t - , j;, m /dr* l ' , )(*l)} Ct.T ■ 

or equivalently 

= j', c (d m j;. m /dc r ”Ko s«(o dt/2<<r j; <m /dO(. i > 

- L k m :,' < K' m * k K'.T ) c 

It remains to compute 

f’i t (d m 0 N \ m /dC m )(O P a (t) dC = - j’, (d m J* +m /d?2 m )(C) P K+1 (0/2(« + 1) dt 
= (l/2(oc+1)) (d m+, J N K +m /dL m+1 )(C) P« +1 tt) dt; 

= ((N + 2oc + m+1)/4(a + D) j* t (d m J N 0t ^ 1 1 tm /dC m )(O P K+ i(0 d£ 

= ((N + 2a + m + 1)/4(a + 1)) . 
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Setting now T tt,m = (N + 2cx + m +1) X^ 1 1,m /8 (a +1) (d m JjJ +m /dt; m )( + 1) and applying once more 
Corollary 11.1 and (11.9), we obtain the result. 


Remark 11.5 : Of course, applying formula (11.10) to the function = 1 , one could think to 
compute the weights p£’™ by 


(11.33) 


a,m a,m 0 

<? 0 .- = <?o. + = 2 


2« 


r (a + 1 ) 2 

l~(2oc + 2) 


- (i/2) z: j= i Qj 01 ' 1 


But, due to the round-off errors, this formula is not so precise as (11.29) and (11.31) for large 
values of N. As a matter of fact, computing pg'™ by using the recursive formulas (11.28), (11.30) 
and (11.32) is also cheaper and easier. 


Remark 11.6 : The computation of the weights from the preceeding formulas seems a little bit 
complicated. However, we intend to work with low values of m ! Tables giving these weights for 
different values of the parameters a, m and N can be found in Appendix C. 

We explicit here the weights when m is equal to 2, since we will use the corresponding 


formula in the sequel : the weights p“' 2 , 1 < j < N, are given by 


(11.34) 


« 2 

pj = 


2 2tt+5 l~(N toe-t2) r(N+0<43) 
N! r(N + 2oc + 5) 


1 


( i - ar 2 ) 2 ) 2 ^^ 2 

the weights Pg' 2 and p“' 2 are given respectively by 


K +2 '(t‘ +2 ) 


(11.35) p 0 K ; 2 = p 0 K ; 2 = 2 


2ou2 r<oe4-1) r(0C43) 


N! 


and 

(11.36) 


p*;? = - p“;! = 2 2k+2 r(a + 2) r(a>3) 


a + 3 r(N + 2 a +5) 

{(a + 2) N 2 4 (oo2)(2a+5) N 4 (<x+3)(2oc+ 3)} 

N! 


r(N + 2a + 5) 

Thanks to the Stirling's formula, it is an easy matter to note that, when N goes to oo, the p*'™ tend 


to 0 as N" 2(Ua) ^ 


-2(2 + a) 


and the p,^ , 0 < k ^ m- 1, tend to 0 as NT 
Let us consider the special (and simpler) case a = 0. Then, the weight p„ is equal to 1 and the 

I R 

Jacobi polynomials simply coincide with the Legendre polynomials. The weights pg' 2 and p°' 2 are 
given by 

(11.37) 
and 


8 


o 0,2 □ o 0,2 = _ 

Po -- P °' + 3(N+1)(N + 2)(N + 3)(N+4) 


(2N 2 4 ION 4 9) 
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(11.38) 


0,2 „ 0,2 
Pi.- D “ Pl.4 = 


8 


(N + l)(Nf2)(N + 3)(N + 4) 


Remark 11.7 : There exists another standard formula to approximate <P(C) p K (C) dt;, namely 
the (left) Gauss-Radau formula 

e.(c) dc - Ej".i^(c“ GB ) e“°” ♦ *<-!) e ! 60 , 
which is not of type (11.10). However, both the Gauss-Radau formula and formula (11.10) are 
special cases of a last formula that we now introduce. 

Let us denote by ( J“ ,p ) n6W the family of Jacobi polynomials which are orthogonal with 
respect to the measure (1 +O a ( 1 -O p dt;, where J“ ,p has degree n. For a fixed integer N, we denote 
by L" p , 1 < j ^ N, the zeros of Jjj ,p . Then, for any integers m ^ 0 and p ^ 0, there exists a 
unique vector of weights p“' m,p , 1 < j < N, and g“'™' p , 0 < k < m -1, and p“'™' p , 0 ^ f < p-1, in 
[R n x IR rrup such that the quadrature formula 


(11.39) 




* l ::', 1) * l ’ i'o (c^/dc'x * i > . 

is exact on P 2N+m+p _i(A). Clearly, the formulas (II. 10) coincide with formulas (11.39) in the 
case m = p, while the (left) Gauss-Radau formula is obtained for m □ 1, p = 0. 


i 
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lll. Variational _ formulation of the. M r.m M . e t problem—for—iim 

bi l g pl a c l a DL 

We are interested in the approximation of the following model problem : Find a function u 
defined on A such that 

u (,v, = f in A , 

(III. 1) u(± 1) bO , 

u'(± I) bO , 

where f is a given distribution on A. It is well-known that, for any f in H” 2 ( A), this problem has 
a unique (variational) solution in Hq(A). However, we need here a more general formulation. 

To that aim, for each parameter cx > - 1, we introduce a family of Sobolev spaces associated 
with the measure p K (C) dC , where the weight g K has been defined in (II. 1). First, we introduce the 
space 

(II 1.2) L 2 ( A) = { v : A-+IR measurable ; v 2 (0 P«(0 dC <+oo } ; 

it is a Hilbert space for the scalar product 

(111.3) (u.v)^}’, utt) v(0 0.(0 dt , 

and we identify it with its dual space; the corresponding norm is noted II• llo.ot.A • ^hen, f° r any 
integer k ^ I, we introduce the Sobolev space 

(111.4) H k (A) = { v € L 2 (A) ; d'v/dt* € L*(A), 0 k } ; 

this space is provided with the norm 

(111.5) ||V|| MA = [|;, Zm <dWV 9.(0 do’' 2 
and with the semi-norm 

dll.6) |v| kotA = [j’, (d k v/dt. k ) 2 p tt (Odt;] ,/2 . 

For any real number s ^ 0 which is not an integer, the space H*(A) is defined by interpolation 
between H^^A) and H^ ]+1 (A), where [s] is the integral part of s, and its norm is denoted by 
||.|| socA • For any nonnegative integer k, we denote by H k 0 (A) the closure in H k (A) of the space 
£)(A) of all functions of class (5°° having a compact support in A; we call H" k (A) its dual space, 
which is spanned by the derivatives of order k of all functions of L 2 (A), and we still use the 
notation (,,.) K for the duality pairing between this space and H k 0 (A). 
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A detailed study of the properties of such spaces can be found in [Gr][BM2], we just recall 
some of them that we need later; assuming that oc satisfies -1 < a < 1, we have 

1) for any k ^ 0, the mapping : v -> vp ft is an isomorphism from H k 0 (A) onto H k K0 (A); 

2) for any integer k ^ 1, the trace mapping : v -> ( v(+ 1), v'(± 1), .... (d k_1 v/dt k_1 )(± 1) ) is 
linear continuous from H k (A) onto IR 2k , and its kernel is exactly the space H k 0 (A); 

3) for any integer k > 1 , the semi-norm |.| kot A is a norm on H k 0 (A), equivalent to ||.|| kct A . 

Next, we consider the bilinear form a K , defined on H 2 (A) x H 2 0 (A) by 

(111.7) a tt (u,v) = u"(U (vp a )"(0 dt; . 

Clearly, for any f in H^ 2 (A), problem <111.1) is equivalent to the following variational one : Find 
u in H 2 0 (A) such that 

(111.8) Vv<eH 2 0 (A), a K (u,v) = (f,v) K . 

In order to study this problem, we need some properties of the form a^. 

The following lemma provides two extensions of the Hardy's inequality that will be useful in 
the sequel. 

Lemma III. 1 : For any real number §, every function <p in £>(A) satisfies 

(111.9) <p ,2 (0 p p (t;) dt; ^(1-2g) <p 2 (0 pp_ 2 (0 dC 

and 

(111.10) \[ y <p ,2 (t;) C 2 pp(OdC >( 1-2^1', <p 2 (Ot 2 pp_ 2 (OdC • 

Proof : The inequality (III.9) follows from 

CK [«p'(t;) pp(C) + (1 -2&) <pU) C Pp.,(0] 2 P.p(0 dC 

< 1 1 , [<P' 2 (0 0p(O + (1 -2&) 2 <p 2 (0 C 2 Pp. 2 (0 + (1 -2&)(«p 2 )'tt) K Pp./O] dC 

^ ( l , U' 2 (0 pp(0 dt; . (1 -2&) l , <p 2 (0 Pp. 2 (U [(1 -2&)C 2 - (1 -K 2 ) + 2t 2 (e-1)] dC 

< (’, «P' 2 (0 Pp(U dt; - (l-2&), <P 2 (0 pp. 2 (0 dC 

Similarly, we also have 

CU J’, WO K PpU) +<p(UPp. 1 (C)] 2 P.p(OdC 

^ j’, [«P' 2 (U t 2 Pp(0 ^ «p 2 (t) P p . 2 (0 + (f 2 )'(U C Pp.,(U] dC 

< <P' 2 (0 C 2 Pp(C)dt; + |’ 1 «P 2 (0pp. 2 (0[1 - (1 -C 2 ) + 2C 2 <&-D] dC 

< J^, <P' 2 (0 t 2 Pp(C) dt; - (1-2&) (’, (P 2 (0 C 2 Pp. 2 (C) dC 
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ProposHion 111.1 : Let cx satisfy -1 < <x < 1. The form a a is continuous on 
H 2 (A) x H 2 0 (A) end elliptic on H 2 0 (A). 

Proof : The continuity of a K is an immeeiiate consequence of the inequality 

K(u,v)| ^ »u"|| 0 a A ti(vp tt )"|| 0 . KiA < liu|| 2ecA l|vpJ 2 _ ttA , 

0 0 

since the mapping : v -* vp R is continuous from H* 0 (A) onto H_ k 0 (A). 

To study the ellipticity, we first note that, for any u in H 2 0 (A), setting w = up K , we have 
a^u.u) = a_ a (w,w) and l|u || 2 a A < c ||w || 2 KA . 

Hence, it is sufficient to prove the ellipticity of a K for a ^ 0. Next, we compute for any a, 
- 1 < a < 1 , 

a K (u,u) = J’, U" 2 (0 p tt (U dK + 2 u”(0u*(0 p;(0 dt > u”(Uu(0 p;(U dC 

= j’, u" 2 (t) P a (U dC - u , 2 (C) p;(0 dC 

~ l-i Pitt) dC - u'tt)utt) p [ ni) (0 dC , 

whence 

(Hl.ll) a„(u,u )={’ 1 u" 2 (U p^COdt; -2j’, u ,2 (U p;(t) dt 4 (t/ 2 ) u 2 (0 pi lv> (Odt. 

We need 

(111.12) p;(C) □ P a _ 2 (0 (-2a) [1 4 (l-2cx) C 2 ] , 

(111.13) p£ v) (0 = p K . 4 (U (-da)d-a) [3 + 6(3-2cx) K 2 + (l-2a)(3-2a) t, 4 ] . 

We must study separately two cases : cx ^ 1/2 and 0 < cx < 1/2. 

1 ) I n the case a ^ 1 / 2 , we have 
(1/2) U 2 (C) pi ,v) (0 dC 

= -2oc(1-oO J^ t u 2 (U P K . 4 (C) [3 4 6(3-2cx) t 2 ♦ (1-2oO(3-2cx) K 4 ] dt 
^ -2a( 1 -a) u 2 (C) P*_ 4 (0 [3 + 6(3-2cx) C 2 ]dC . 

Using (III.9) and (111. 1 0) for 0 = cx- 2 , together with a density argument, we derive 
( 1 / 2 ) J’, u 2 (t) pi ,V) (U dC 

^ - [2cx( 1 -cx)/(5-2cx)] u' 2 (0 p K _ 2 (U [3 + 6(3-2cx) K 2 ] dt; 


That yields 
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a tt (u,u) s* u" 2 (U q k (U dt; + [2oc/(5-2oc)] u' 2 (0 p w _ 2 (t;) 

[2(5-2a) + 2(5-2a)( 1 -2a) t 2 - 3( 1 -a) - 6( 1 -a)(3-2a) t 2 ] dt, 

^ f-i u" 2 (t) 0«(C) dC + [2a/(5-2a)] j], u' 2 (U p a _ 2 (0 

[(7-a) - 2(4-3a + 2a 2 ) t 2 ] d£ 

Noting that, for ang C in A and a > 1/2, 

(7-a) - 2(4-3a + 2a 2 ) C 2 > (7-a) - 2(4-3a + 2a 2 ) n - (1-5a + 4a 2 ) 

= - (1-a)(1-4a) > 0 , 

we finally obtain 

o a (u,u) ^|u| 2 kA • 

Since |.| 2 K A is equivalent to the norm I!• ll 2 « a on H 2 0 ( A), that gives the result. 

2) In the case 0 a ^ 1/2, using again (II 1.9) and (111.10) for £ = a-2 and a density argument 
yields 

a K (u,u) ^ u" 2 (C) P K (U dt; + 4a(5-2a) f ^ u 2 (0 p k _ 4 (U [ 1 + (1 -2a) t; 2 ] dt; 

-2a( 1 -a) u 2 (0 p K _ 4 (C) [3 + 6(3-2a) t, 2 + (1-2a)(3-2a) K 4 ] dt; , 

that is to say 

a ot (u,u) > u ,,2 (0 P K (0 dt; + 2a f A u 2 (t) P a _ 4 (0 

[(7-a) - 2(4-3a + 2a 2 ) K 2 - (1 -a)( 1 -2a)(3-2a) C; 4 ] dC . 
We note that, for any t; in A and a ^ 1/2, 

(7-a) - 2(4-3a 1 2a 2 ) t; 2 - (1-a)( 1-2a)(3-2a) t; 4 

^ - (1 -a)( 1 -4a) - (1 -a)( 1 -2a)(3-2a) = - 4( 1-a)( 1-3a+a 2 ) , 

which proves the ellipticity for a > (3-V5)/2. 

In the case a < (3-%/5)/2, we fix a constant X, 0 < X < 1, and we use (111.9) for £ = a 
then for £ = a-2 , again with a density argument, to estimate the first term. We obtain 
a K (u,u) >(1-X) u" 2 (0 p K (Odt + X( 1 -2a)(5-2a) }’, u 2 (U p tt . 4 (0 dt; 

- 8a(l-a)(1-3a.a 2 ) u 2 (0 P tt . 4 (0 dt . 

In order to choose 

i 

X = 8a( 1 -a)( 1 -3a+a 2 )/( 1 -2a)(5-2a) , 

we just have to check that, if a satisfies 0 ^ a < (Z-V5)/2, the right-hand member of the 
previous line is < 1. This last condition is equivalent to 
8a( 1 -a)( 1 - 3a +a 2 ) < (1 - 2a)(5-2a) , 

i.e. to the positiveness of the polynomial 
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P(a) = 8a 4 - 32a 3 + 36a 2 - 20a + 5 . 

Computing P"(a) = 24(4a 2 - 8a + 3) = 24( 1 -2a)(3-2a), we see that, when a goes from 0 to 
1 /2, P' increases from -20 to -8 and P decreases from 5 to 1/2, hence it is > 1 /2. That ends the 
proof. 

Remark 111.1 : In the particular case a = -1/2 of the Chebgshev weight, the properties of a K 
have been proved in [M, Lemma V.1 ] by a slightly different argument. 


An immediate consequence of Proposition III. 1 is the 
Theorem 111.1 : Let a satisfy -1 <a<1. For any f in H~ 2 ( A), problem (II 1.1) has a 
unique solution u in H 2 0 (A). Moreover, it satisfies 

Cii.m) MI 2 ...AsUlf|l„ ;J|A1 

We are also in a position to propose a first approximation of problem (111.1). Let N be a 
given integer. The discrete problem is the following one : Find a polynomial u N in 
P n (A) Pi H 2 0 (A) such that 

(111.15) V v N € P n (A) fl H 2 0 (A), a K (u M , v N ) = (f,v N ) a . 

Theorem 111.2 : Let a satisfy -1 <a< I. For any f in H^ 2 ( A), problem (111.15) has a 
unique solution u N in P N (A) fl H 2 0 (A). This solution converges to u when N tends to +co. 
Moreover, if the solution u of problem (lll.l) belongs to H*(A) for a real number as* 1 , the 
following error estimate is satisfied 

(111.16) llu-u N || 2 w A <cN 2 " T ||u|| (r K A . 


Proof: Proposition 111.1 Implies at once that problem 
also yields that, for any w N in P N (A) fl H 2 0 (A), 

II u n- w nII|,oc.a ^ ca «< u N- w N - u n-w n ) =ca,(u-w 


(111.15) is well-posed; due to (111.8), it 

N > U N -W n) ^ C II U_W nH2,ci,A II U N W N ^2,«,A ’ 


whence 


(111.17) llu-u N || 2(tA < c inf H u-w nII 2 .«.a • 

w N € P n (A) 0 H 2 0 (A) 

Choosing w N = n^' 2 u, where n^' 2 is defined in Appendix A, we deduce (lll.l 6) from Theorem A.6. 
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The convergence is then derived bg a standard density argument. 


Remark 111.2 : By an interpolation argument, one can easily prove that, whenever the data f is 
in the space H®(A) for a real number p ^ 0, the solution u of (lll.l) belongs to H£ +4 (A). 
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IV. Two collocation methods for the Pi rich let problem _fa r ..t he 

bilaplacian. 


The aim of this section is to analyse and compare two different collocation methods to 
approximate problem (1(1.1). The advantage of these methods is that the corresponding mass 
matrix is diagonal, which is of importance for some algorithmic reasons : the explicit resolution of 
time-dependent problems or the design of simple preconditioners for instance. It is well-known 
that, if one wants a spectral collocation algorithm to be accurate, the collocation points must be 
chosen as the nodes of a quadrature formula. Hence, these methods are related to the quadrature 
formula (11.10) with respectively m = 1 and m = 2. 

Let us introduce some notation. For any integer m 0, we consider the bilinear form 
(,,.) Km defined on 0 m_1 (A) x C m ' , (A) by 

(tv. i> («,*).„ = e",, »(cr m ) *<er m > or * £r:» (A«)/deKi» ■ 


(IV.2) 


m 


Now, we present the discrete problems. For a function f continuous on A, they are the 
following ones : Find a polynomial u N in P Nt3 (A) such that 
u‘ lv, (t;“ +m ) = f(C; +m ) , l<j<N , 

U N (±l)rO , 

U’(±1) = 0 , 

where m is equal either to 1 or 2. In both problems, the number of equations is equal to the 
dimension of P N+3 (A). 


We begin by studying problem (IV.2) 2 , which will turn out to be the easier one. We have 
the 

Proposition IV. 1 : Let a satisfy -1 < a < 1. Problem (IV.2) 2 is equivalent to the 

r\ 

following variational one: Find u N in P N+3 (A) f) H k0 (A) such that 
(IV.3) Vv N €P N+3 (A)nH K 2 0 (A), a a (u N ,v N )u(f,v N ) a2 . 

Proof : By multiplying the first equation in (IV.2) 2 by the (1 -L 2 ) 2 Q* +2 , 1 < j < N, which form 
a basis of P N+3 (A) H 0 ( A), we see that (IY.2) 2 is equivalent to find u N in P N+3 (A) fl 0 ( A) 
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such that 

^ V N € ^N+3^^ H*<j(A) > ^ U N > V N^«,2 = ^> V N^oc,2 

Next, since the quadrature formula (11.10) for m = 2 is exact on P 2N+3 (A), we observe that 
(C> . v N ) tt2 * s ec l ua1 (° (u‘ IV) - v n)* : integrating by parts gives the result. 


As in Section III, we obtain the 

Theorem IV. 1 : Let oc satisfy -1 < a < 1. For any f in C°(A), problem (IV.2) 2 has a 
unique solution u N in P N+3 (A) fl H* 0 (A). 

We want now to estimate the error between u and u N . 

Theorem IV.2 : Let oc satisfy -1 < ot < 1. If the solution u of problem (lll.l) belongs to 
H*(A) for a real number a ^ 1, and if the data f is such that the function (I -t 2 ) 3/2 f belongs 
to a space H^(A) for a real number q > 1/2 , the following error estimate between the solutions 
of problems ( 111.1) and (IV.2) 2 is satisfied 
(IV.4) ||u-u N || 2 , A < c ( N 2 -' ||u||, „ A * n' /2 -> ||( I -t 2 ) 3 ' 2 f|l, x . A ) . 


Proof : Using Proposition 111.1 together with (111,8) and (IV.3) gives, for any w N in 

K" W jf,«.A < Ca ^ U N- W N . U N- W N> = Ca ^ U " W N • U N* W n) - ( f ' U N- W N) R f ( f - U N- W N)oc,2 • 

whence 


(IV.5) 


|U-U N |( 2 A < C( inf 


w„«P N . 3 (A)flH. 2 „(A) 




Vo v 


sup 


v»« n h 2 , 0 <a) 


(f,v N ) K -(f,v N ) cc , 2 

II V N ll2,<x.,A 


) 


We choose w N = n^jU; then, we deduce the result from Theorem A.6 and Lemma B.2 of the 


appendices. 


Remark IV. 1 : The smoothness assumption we make on f is very weak, since we do not require 
that f is continuous on A, but only on A. 
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Now, we study problem (IV.2), . Here, we must define the following bilinear form on 
0 4 (A) x 0°(A) 

(IV.6) a K ,N (u - v) = (u<IV) - v) «,t • 

This form no longer coincides with the a ot (.,.) on P N+3 (A) D H* 0 (A). Nevertheless, we have 
Proposition IV.2 : Let a satisfy -1 <a< I. Problem (IV.2), is equivalent to the 
following variational one: Find u N in P N+3 (A) f) H* 0 (A) such that 
(IV.7) V V N € ^N+3^^ H*o(A), 1 V N^ = ^ ,V n)k,1 

Proof: We obtain (IV.7) simply by multiplying the first equation in (IV.2), by the (1 -t 2 ) 2 Q“ +1 
, 1 ^ j N, which also form a basis of P N+3 (A) f] H* 0 (A). 


Studying the form a KN requires the following lemma. 

Lemma IV.2 : For any real number oc > -1, the following inequalities hold 
(IY.8) c 1 II^N+t llo,«,A ^ ^N-1 '^N+3^,1 ^ C II^N + 1 llo,«,A ' 


Proof : Let us introduce the following notation : given two quantities X(N) and y(N) depending on 
N, X(N) ~ y(N) means that there exists a positive constant c independent of N such that 


c" 1 X(N) < y(N) c X(N). Using the induction formula (11.5) with n = N + 2, we compute 

( J N-1 1 J N + 3^oi,1 = 

(2N + 2cx + 5)(N + cx + 3) (£ J*_, , ^n +2 )«.i " (N+oc + 2)(N+oc +3) (J*., , J* t1 ) K l 

(N + 3)(N + 2a + 3) 

but, since (J^_, , JU +1 ) K , coincides with (, , J|J +1 ) K which is equal to 0, we simply have 

("Vi + '^N+2^,1 

The same formula (11.5) applied with n = N-1 to compute t; J^_, gives 

(^N-1 ' ^N*3^<x,1 • ’ ’ ^N + 2^*.1 

whertee, by (11.6), 

^N-1 >^N + 3^«,1 ^N + 1 )J N + 2^«,1 

Since the nodes of formula (II. 10) for m = 1 are the zeros of , we deduce 

^N-1 > ^N + 3^«,1 ~ N 1 Po,± + 1 ^N + 2^ ^ • 

But, due to formula (11.8), we have 
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2(a+1) J*'(±1) = n(n + 2a +1) J*(± I) 

and, from (11.3) and the Stirling's formula, we derive that, as n goes to + oo, J*(± 1) behaves as n a 
up to a multiplicative constant. Hence, using (11.20) and once more the Stirling's formula, we 
obtain 

(Jn_, . Oo ~ N '’ N ' 2 ' 2 * n2+ * N* ~ N - 1 . 

On the other hand, it follows from (11.4) and the Stirling's formula that ||Jn + 1 IIo,«,a ~ N_1 - 
whence the lemma. 


Corollaru IV. 1 : For any real number oc > - 1 .the following inequalities hold 


(IV.9) 


N 4 m„\, \\l.. A < w„\r ■ 0 n \ 3 ).., <«N 4 k . , ii„ 2 ,. a 


Proof: Derivating four times formula (11.6) with n □ N + 2, we see that 
,oc <iv) _ (2N + 2a + 5)(N*2<x + 3) |0t (m) _ N+oc+2 (iv) 

N+3 ” N + oc + 3 N+2 + N + 2oc + 2 N + 1 

Repeating this argument three times more, we obtain 

°N+3 - a n °n-i + r N > 

where X N is a real coefficient, \ N ~ N 4 , and r N is a polynomial with degree < N-2. Then, (IV.9) is 
a simple consequence of (IV.8). 


We can now prove the 

Proposition IV.3 : Let <x satisfy -1 < oc < 1. The form a aN satisfies the following 
properties of continuity 

(IV. 10) V u N € P n+3 (A), V v N € P n+ 3 (A) D H 2 0 (A), |a KN (u N ,v n )|<c||u n || 2aA ||v n || 2otA 
and of ellipticity 

(IV.ll) V U N € fl H 2 0 (A) , a «,N^ U N > U N^ ^ C II U N il2 ,«.,A 


Proof : Let us write any polynomials u N and v N of P N+3 (A) on the form 

L N + 3 „ n - , N + 3 ^p * 

n = 0 U an< ^ V N = *-nrO V 

Since the quadrature formula (II. 1 0) with m = I is exact on P 2 N+ 1 ^^^• we derive from (III.7) 
and (IV.6) 

a «,h/ U N ’ V n) " a «^ U N ’ V N^ = ^ ^ ^N + 3 ,v) n + 3^(x,1 
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This formula, together with Proposition 111.1 and the positivity of (J^ 17 * , ^+3^*1 P r oved in 
Corollary IV. 1, yields the property of ellipticity (IV. 1 1). On the other hand, in order to prove 
(IV. 10), it suffices to check that for any v N in P N+3 (A), 

(IV.12) (v ) (J N+3 * .^ c l V J2,«.A ■ 

Writing now v^' = Lp : + o z" J“ , we have 

lv N l 2 2 , tt , A = L n N :o (zV||jX.a ; 

using (11.6), we compute 

Jjf + 3 = (2N42oc + 3)(2N + 2a + 5) 


(N+oc + 2)(N+o(43) 


'^N+l ' 


(N< 2a + 2)(N + 2a + 3) 
where the degree of q N is ^ N; hence, comparing the two expansions of v^' yields 


z N+1 = (2N + 2oc + 3)(2N + 2<x + 5) 


(N + ix + 2)(N + a + 3) 
(N + 2a + 2)(N + 2cx + 3) 


'N+3 


so that 


IvJ 2 z . k . a ^on‘’(v n -) 2 ||j“ + 1 || 0 V a • 

This inequality, together with Corollary IV. 1, implies (IV. 12), hence the proposition. 


As previously, we derive that the discrete problem is well-posed. 

Theorem IV.3 : Let <x satisfy -1 < oc < 1. For any f in 6°(A), problem (IV.2), has a 
unique solution u N in P N+3 (A) D H^ 0 (A). 


The error between u and u N is given in the 

Theorem IV. 4 : Let a satisfy -1 < a < 1. If the solution u of problem (111.1) belongs to 
H"(A) for a real number a ^ 1, and if the data f is such that the function (1-C 2 )f belongstoa 
space H£(A) for a real number p > 1/2, the following error estimate between the solutions of 
problems ( III. I) and (IV.2), is satisfied 

(IV, I 3) ||u-u„ Il 2 „ A < c ( N ! - Hull, „ A . N ,/2 -' ||( I -t 2 ) f||,„ A ) 


Proof : Using Proposition IV.3 together with (III.8) and (IV.7) gives, for any w N in 

P N.3< A >nH 2 0<A), 

HU H - W N III...A < Ca «,N< U »- W N ■ U M- W »> ='«>-»« • " N - W n) - < f ’ U lT W N>« * 

so that 




- 27 - 


(IV. M) II U_U N H2.cc, A ^ C ^ inf II u— W N l^ ot.A 

W N ^ ^N+3^A) f~l H m (A) 

(f,v N ) #t -(f,v N ) 1 

+ sup -—-— ) 

v N €P N+ 3 (A)nH^ 0 (A) II v nII 2 ,oc.a 

We choose w N = n^u; then, we deduce the result from Theorem A.6 and Lemma B.2 of the 
appendices. 


Remark IV.2 : Here again, the continuity of the data f on A is not necessary for the discrete 
solution to converge to the exact one. However, the assumption of Theorem IV.4 is stronger that the 
assumption of Theorem IV.2 for the same order of accuracy : indeed, taking f(0 □ (1 we 

know that (1 -t; 2 ) 3/2 f belongs to H®(A) if and only if g Is < § + (4 tcx)/2, while (1 -C 2 ) f belongs 
to H 9 (A) if and only if p is < § + (3+<x)/2. 


Remark IV.3 : By applying a standard duality method, one could obtain an error bound in the 
norm IUI 0(xA . Precisely, under the assumptions of Theorem IV.2 or IV.4, it would be possible to 
prove the following estimate between the solution u of problem (lll.l) and the solution u N of 
problem (IV.2) m 

(IV. 1 5) l|u-u N || 0 ct A < c ( hr* Hull, kA + N 1/2 " 9 1|( 1 -t 2 ) (m+,,/2 f II ? , k . a ) ■ 

However, the term N 1/2 " 9 1|( 1 _£ 2 ) <m+1,/2 f|| is not improved; hence, this last estimate is not of 
great interest, since this term is the worst one (indeed, the fact that (I -C 2 ) m f belongs to H 9 (A) 
implies that f is in H 9 " m (A), hence that u belongs to H 9+4 ~ m (A), and the other term in (IV,4) or 
(IV. 1 3) is N m-2 "° ||u|| p+4 _ mK A , which is smaller). 


Remark IV.4 : Several other collocation methods seem natural to approximate problem (lll.l). 
Let us consider two of them and prove that they are not so good as problems (IV.2), and (IV.2) 2 . 

First, an immediate analogue of these problems consists of using the Gauss points, which 
gives: Find a polynomial u N in P N+3 (A) such that 


C’ttp = rap 
(IV.16) | u N (± 1) = 0 , 

uj/± 1) = 0 . 

Of course, setting now 

a wN (u,v) = (u (lv) ,v) tt0 . 


1 < j <N 
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we see that this problem is equivalent to the following variational one : Find u N in 
P N+3 (A) H H 2 q(A) such that 

Vv ^U A >nH>), ,v N ) = (f,v N ) K0 • 

However, choosing v □ (1 -t; 2 ) 2 " , we observe that 

o aN (v,v) c N 6 and ||v|| 2 cc A > c' N 7 , 

which shows that the constant of ellipticity of the form a^ N is not bounded from below 
independently of N. Consequently that the approximation of u by u N is not optimal. 

Another method consists of searching u N as a polynomial of lower degree, which is achieved 
in the following : Find a polynomial u N in P N+1 (A) such that 
uJ v) (Cj t * m ) = f(Cj a+m ) , 2 < j < N-1 , 

(IV.17) u n (±1)dO , 
u^(± 1) = 0 , 


where m can be equal to 0, 1 or 2. This problem admits the variational formulation : Find u N in 
P N+ ,( A ) fl H^ 0 ( A ) such that 


^ v n € p n+i^ a ^ v n^i ) - v n^n ) - 0- ( u n * > - ^ ,v N ) Km 

But, in this case, the norm of the form : (u,v) -> (u^’.v)^ on the discrete spaces is no longer 


bounded independently of N, which prevents the discrete problem to be well-posed (indeed, in the 


simplest case a □ 0, denoting by k the integral part of (N-1 )/2, we set 
u K «) = ft 4,.,:XOdt and v N (t> = (1 -C 2 K(C7'> 2 -C 2 ) 

it is easy to check that 

H » v n^o I ^ c N and II u n H 2 ,o,a ^ c ^ * 

so that the norm is ^ c N 1/2 ). 
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V. A collocation method for the periodic nonperiodic Dirichlet 
problem fo r t h e, t? 1] gp|gc1an> 


We complete the analysis of the best collocation method described in Section IV (collocation 
at the N nodes of the generalized Gauss quadrature formula for m = 2) by applying it to a 
bidimensional Dirichlet problem for the bilaplacian, when the boundary conditions are periodic in 
the first direction and homogeneous in the second one. More precisely, let 0 denote the domain 
0 x A , where 0 and A stand respectively for the intervals ]-n ,n[ and ]-1,1 [. The generic point 
in 0 is denoted by x = (x ,y). 

For a given distribution f on 0, the problem we want to approximate is the following one : 
Find a function u defined on 0 such that 
I A 2 u = f in 0 , 


(V. 1) | u(-n,y) = u( + n,y) and u(x,± 1) = 0 ,x = (x,y)€Q 

(9u/dx)(-n ,y) = Ou/dx)( + n ,y) and (5u/0y)(x,±t) = 0 , x = (x,y)€0 

The boundary conditions are periodic in the x-direction and homogeneous Dirichlet in the 
y-direction. 


Example V. i : An example of problem (V.l) is given by the stationary Stokes equations 
governing the flow of a viscous incompressible fluid between two parallel planes y □ ± ), when the 
body forces are parallel to the plane spanned by the x and y-directions and depends only on the two 
coordinates x and y : 


(V.2) 


- v Au + grad p = f 

div u = 0 


Here, u represents the velocity of the fluid, p is the pressure and the kinematic viscosity v is a 
fixed parameter > 0. A model of this problem is obtained by reducing it to the previously defined 
domain 0 and providing it with the boundary conditions 
(V.3) u(-n,y) = u( + n,y) and u(x,±1) = 0 ,x=(x,y)€0 . 

Then, following the techniques of [BMM, Prop. 11.2], we can prove that, for any u in H^CQ), 
a^l, there exists a unique stream function u in H ff+, (Q) such that 
(V.4) u = curl u in 0 
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This function u is a solution of problem (V. 1) with f = curl f. Similar results can be obtained in 
weighted Sobolev spaces [MM], 


The analysis of problem (V. 1) requires some notation. First, if X is a Banach space, for any 
real number s > 0, we denote by H’(0,X) the closure in H*(0,X) of the space of all functions of 
class 0°° on 0 which are periodic with a 2n-period, and by H^ s (0,X') its dual space. Let v* , 
F <e TL , stand for the Fourier coefficients of any function v in L 2 (0 ,X), we recall that the mapping : 

is a norm on H*(0,X), equivalent to that induced by H*(0,X). Then, defining the Fourier 
coefficients ?*, £ € TL, of any distribution f in H^ s (0,X'), we see that the norm on H" s (0,X') is 


given by 


HjS(0,X*) 


= <L, 


e TL 


a t? 2 r’ 


if ii, 2 .) 


. 1/2 


Next, for each parameter a > -1, for any real number s ^ 0, we introduce the space 
(V.5) H* # (0) = H’(0,L 2 (A)) D L 2 (0,H’(A)) , 

and provide it with the norm 


(V.6) 


Ml,.,.. n = (£, e! << 2 ’ llAl„ 2 .„. A uiP'H, 2 , a 


)) 


1/2 


For any integer k js 0, we also consider the space 


(V.7) 


h^ 0 (q) = h:( 0 ,l 2 (a)) n l 2 ( 0 , h; 0 (a)) 


a,O v 


and we denote its dual space by H‘ k # (0) 


Clearly, if the distribution f is given in H” 2 # (Q), problem (V.l) is equivalent to the 
following variational one : Find u in H 2 #0 (Q) such that 
(V. 8) Vv€H a,„o (Q) ' J-rr (A 2 u,v) a dx = (f,v) K dx . 

Finally, taking the Fourier transform of this equation gives the equivalent formulation : Find u 
in H 2 #0 (Q) such that, for any integer £ , 

(V.9J V z € H 2 0 (A), <fV- 2S> 2 u*" + u ,(IV> , z) K = (f* , z) K . 


0 0 

Let us define the bilinear form on H R 0 ( A) x H a 0 (A) 

(V. 1 0) b*(w,z) = \[ y { F 4 w(y) z(y) p K (y) *■ 2 ? 2 w'(y) (zp tt )'(y) + w"(y) (zp a )"(y) } dy 
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: Let cx satisfy -1 < cx < 1. For any integer £ , the form b* satisfies the 


following properties of continuity 


(V. 1 1) »w(H;,(A),»2<H',(A), 

|b'(w. 2 )| < c (f 4 ||w|| 2 „ A * ||w|| 2 , A ) m ( r II4„ 2 .,, A * l|z|| 2 . A ) 

and of ellipticity 

(V.12) Vw€Hj 0 (A), b*(w,w) > c ( £ 4 llw||o i(lA + || w||| ttA ) 


2 sW2 


Proof: Integrating bg parts, we have for any w and z in H 2 0 ( A), 

b*(w,z) □ \[ y {£ 4 w(y) z(y) p a (y) - 2t 2 w(y) (zp K )"(y) + w"(y) (zp tt )"(y)} dy , 
so that 

|b'(w,z)| < £ 4 l|w|| 0 K A ||z|| 0a a + 21 2 ||w|| 0r A l|zpJ| 2 . a A + l|w|| 2 a0 l|zpJI 2<0tiA . 

Using the continuity of the mapping : v -> vp K from H 2 0 (A) onto H 2 k 0 (A), we obtain (V. 1 1). 
Moreover, for z = w, we write 

b*(w,z)n j’, {f 4 w 2 (y) p^Cy) + 2£ 2 w’(y) (wp ct )'(y) +w"(y) (wp ec )"(y) } dy , 
and (V. 12) follows from Proposition II 1.1 and the ellipticity [BM2, Lemma 111.4] of the form 
(w,z) -> w’(y) (zp K )'(y) dy on hJqCA). 

Theorem V. 1 : Let cx satisfy -1 < oc < 1. For any f in H" 2 # (Q) , problem (V.1) has a 
unique solution u in H 2 #0 (Q). Moreover, it satisfies 
(Y.I3) Il u ll2,ct,#,0 ^ c ^Hhj2 # (o) 


Proof : Any distribution f in H’ 2 ^(0) can be written f = f # + f 0 , where f, belongs to 
H; 2 (0,L 2 (A)) and f 0 belongs to L 2 (0,H^ 2 (A)), with 

H f llH 5> 2 # (0) < i f # Hh>2(0,l£(A)) + H f 0 Hl2(©,H 5 2(a» ■ 

Due to Proposition V. I, for any integer £ , equation (V.9) has a unique solution w* satisfying 

, « 4 llw'll 2 . A * ||w’|l!„ A < c ((1 <( 2 )- 2 ii?:i 2 ,. a * II?J|I„V<a> ) .. 

The function u, the Fourier coefficients of which are the w* , £ € 7L, is the unique solution of (V. 1) 
in H 2 n (Q) and satisfies (Y. 1 3). 


In order to approximate problem (V.1), we fix an integer L ^ 1. Then, for any Banach space 
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X, we consider the space $ L (0,X) of trigonometric series of order L on 0, i.e., 

(V.M) S L (0,X) = { v L □ v* exp(tfx), v* € X } 

Next, we introduce the nodes 

(V. 1 5) x k = 2kn/(2L+1) , -L ^ k < L . 

Let {,£ denote the interpolation operator at the nodes x k , -L k < L : for any function v in 
0°(0,X), v'v belongs to S K (0,X) and satisfies 
(V. 1 6) (t L # v)(x k ) = v(x k ) , -L < k < L . 


Now, let 6 denote the discretization parameter (L,N), where L is > 1 and N is ^ 3. We define 
the space of discrete solutions as 
(V. 1 7) X 6 = S l (0,P n+3 (A)) • 

The discrete problem is the following one : Find u s in X 8 such that 

-L < k < L, 1 < j < N , 


(V, 1 8) 


A 2 u 5 (x k , L“ +2 ) = f(x k , tf 2 ) 


u.CXfc , ±1) = 0 


-L < k < L 


(du 6 /0n)(x k , ± 1) = 0 , -L < k < L . 

The number of equations is equal to (2L+ 1 )(N + 4), which is the dimension of X s . 


To analyse this problem, we define the discrete bilinear form on 0 S (A) x C’(A) by 
(V. 1 9) b* N (w,z) = (w,z) k2 - 2 ? 2 (w",z) K(2 + ( w<iv, ,z) k>2 • 

We have the 

Proposition V.2 : Let a satisfy -1 < oc < 1 . Problem (V.18) is equivalent to the following 
variational one: Find u 6 in $ L (0,P Nt3 (A) fl H 2 0 (A)) such that, for any integer £ between -L 
and L , 

(V.20) Vz,«PjA)nH>), b' N (u ( ',z N ) = (l.'f',Z H )„ 2 . 


Proof: Since the interpolation operator v* is equal to the identity on 3 L (0,P N+3 (A)), problem 
(V. 1 8) is clearly equivalent to the following one : Find u 5 in S L (0 ,P Nt3 (A)) such that, for any 
integer £ between-L and L, 

G'(t“ +2 ) *2f 2 u^ttj" 2 ) + Uj <lv) (t;* +2 ) = (^f'KCj** 2 ) , 1 < j < N , 

o;<± i > = o , 

(0Gg/0n)(± 1) = 0 


(V.21) 
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Multiplying the first equation in (V.21) by the polynomials (1 -L 2 ) 2 Q* +2 , 1 < j N, we obtain 
(V.20). 


We must now study the properties of the form b* N . We need some lemmas. 
Lemma V. 1 : For any real number o< > — 1 , the following inequalities hold 


(V.22) 


114 


N + 2 


I 2 

"o.d.A 




N+4'k,2 


J* 

; N+2 


|2 

l0,ct,A 


Proof : We follow here the same lines as in the proof of Lemma IV.2 and the notation X(N) ~ y(N) 
again means that the quantity X(N)/y(N) is bounded by positive constants independent of N. Due to 
the definition (IV. I), we have 


(J 


a+2 


> J N+4^«,2 


= 2 J* +2 ( + 1 )J“ +4 ( +1 ) P«; t 2 + 2 (J* +2 ( . 1 )4 + ;( .Of u£ +2 '( +1 )4 +4 ( f i)) p*; 2 . 
From (11.35), (11.36) and Lemma 11.7, we deduce 


(^N > 9-2 


N + 4 / ot,2 


2 k+ 3 r(a+o T(a+3) 
a +3 


N! 


f(N + 2oc + 5) 


4 +2 M)4 +4 ui) 


{ [(a+2) N 2 + (a.2)(2a + 5) N ♦ (a+3)(2a + 3)] 

- (a +1 )(<x + 3) [(N + 4)((N + 2a+5)/2(ocf 1)) + (N(N + 2a + 5)/2(a + 3))] } 


whence 

= - 2 ! "* 5 r<«. I) r(a + 3) r(N<2 N J; 5) J n " 2 M>J n \ 4 (*I)(2N*2o<*7) . 

Using (11.3) and the Stirling formula, we see that - (Jjj +2 > 4+4^2 ~ N" 1 • Formula (II.4), 
together with the Stirling's formula, also yields that II^n +2 Ho «. a ~ N” 1 . whence the result. 


Corollaru V. t : For any real number a > - 1 , the following inequalities hold 
(V.23) c N |U N+2 1| 0A A < - (U N+3 , U n+3 ) k2 < c N l|U N+2 IIo.«.,a 

Prpof : Using twice formula (11.9) implies 

(4+3 - 4+3^ct,2 ~~ N2 (4+? > 4 + 3^.2 

Next, we use the induction formula (11.6) in order to replace U^ 2 by a combination of C4 +2 and 
4-? ■ Not Og that (Jjj* 2 , Un+ 3^2 = (4-i - ^n+3^ ec l ua1 to we obtain 
( J N+3 > 4+3^«.2 ~ (*^4 +2 > ^N + 3^,2 
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We replace t;j£ +3 by a combination of J* +4 and J* +2 , and observe that (J* +2 , J* +2 ) k>2 is equal 
to 0. That gives 


^N + 3 1 vJ N + 3^cc.2 ~ N ' J N + 4^,2 

so that the corollary follows from Lemma V. 1. 


Lemma Y.2 : For any real number oc > - 1 , the following equality holds 
(V.24) ((1 V)J“ " , (1 -L 2 )J n k ") k>0 = 2 (oc + 1) (N-1 )N(2N + 2oc + 1) ||J* || 2 k A 

Proof: Using successively (11.8), (11.5) and (11.6), we write 

(1-C; 2 ) J n k "(C“) = 2 (oc+ 1) tj a J N a 'a“) = 2 (oc + 1) ttJjjmj 1 ) 

= 2 ((oc + 1) /(2N + 2<x + 1)) [((N+ 1 )(N + 2oc+ 1 )/(N+c<4 1)) J«;,(t;p + (N+oc) Jj^O;*)] 
= 2 ((cx + 1)(N.oc)/(N.2cx)) J*;,(t;p 
That gives 

((1 -C 2 )J n k " , (1 V)J N ‘ ") Ki0 = 4 <(oc + 1) (N+cx)/(N + 2cx)) 2 || 2 kA . 

The quantity HJjp, IIq K A can be computed from the Gauss-Lobatto quadrature formula : 

ll J J;Jo,«.A = 2 , js; 1 ( + i) 2 e; : ’ , 

so that, due to (11.20), (11.9), (11.3) and (11.4), 

K;iCa=((N + 1)(N>2«,2)(2N + 2« + 3)/2( ( x + 1)) ||J* +1 || 0 2 kA 

= ((N+ 1 )(N + 2)(N + 2<x+2) 2 (2N + 2oc+5)/2(o( + 1 )(N+cx + 2) 2 ) ||jJ +2 || 2 >otA , 
Using this formula with N+1 replaced by N-1 gives the lemma. 


Lemma Y.3 : For any real number a > - 1 , the following equality holds for any real numbers 
X and y 


(V.25) 


Vex 2 


•’n-i 


■ 2> 


J* 

; N + 1 


'0,0c, A 


) 


^((1 C)(XJ N _ 1 + yJ^' 1 ),(l-t; 2 )(X J^" ( + y Jjj iWo 

< C N (X l|vJ N _ t llo.ec.A + y 2 H J N + 1 IIqa.A ) 


Proof : Thanks to (11.8), (11.6) and (11.5), we have 

(l-c; 2 )^;,^) = 2 (cx+1) - (N+1 )(N + 2oc + 2) J N * +1 U*) 

= 2 (oc + 1) ((N+<x)(N + oc + 1)/(N + 2oc)(N + 2oc + 1)) 

+ ((N+cx)(N + cx+ 1 )(N + 2cx 1 2)/(N + 2oc+ 1)) J^_,(t;p 



- 35 - 


= ((N+oO(N+<x+1)/(N + 2oc)(N+2oc+1)) (1-C* 2 ) 

+ ((N+a)(N+a + 1)(2N + 2a4l)/(N + 2a+D) Jj./tp 

That implies 

((i-c 2 )(xj *: 1 + y j N a ; 1 ) ,d-c 2 )(x j*:, + yj 

= ||[X + y ((N+a)(N+a + 1)/(N+2a)(N+2a+1))] (1 -C 2 )^n-i 

+ y((N + oO(N + oc+1)(2N + 2oc. 1)/(N+2oc+ 1)) ■ 

By using twice formula (11.8) , we obtain at once 

ii<i-t 2 )Ciiio.«,A =I-i 'JN-,(U(P (t+2 j N K : 1 ) ,, a)dc 

= (N-2)(N-1)(N42oc)(N42<x+1)||Jn_, \\ 2 0aA . 

Noting that (1 -C 2 )^!', + (N-l)(N-2) Jj_, is of degree < N-2, we see that 
)', + (N-lKN-2) Jj./O] ^.,(0 P a (OdC = ° , 

hence 

{', P«(C)dC = -(N-2)(N-1)|M*. 1 || 0 2 aA . 

Finally, setting 

y* = y((N+oc)(N+o<*1)/(N4 2oc)/(N + 2oc+1)) , 

we have 

= (X + y*) 2 (N-2)(N- 1 )(N + 2oc)(N + 2o(41) 

- 2 y* (X + y*) (N-2)(N-1)(N42a)(2N + 2oc4l) + y* 2 (N + 2a) 2 (2N + 2a41) 2 
= (N + 2a) { X 2 (N-2)(N-l)(N + 20C4l) - 2 Xy* (N-2)(N- 1 )N 

+ Y* 2 [N 3 + 2(7cx + 6)N 2 + 2(<x + 1)(1 Ooc- 1 )N > 2(2oc + 1)(2oc- 1)(oc + 1)] } . 

Using the inequality 2 [Xy*| < X 2 + y* 2 , we deduce that the left-hand side is less than 
c N 4 (X 2 + y* 2 ), and greater than 

(N + 2oc) { X 2 (N-2)(N-1)(20C4 1) 

. y* 2 [( Mo(4 15)N 2 4 2(10a 2 49a-2)N . 2(2oo1 )(2a-1 )(oc + 1)]} . 

It is easy to check that this last polynomial can also be written 
(1 40(4 15)N(N-2) 4 2(0(4 1)[(10oc+ 13)N 4 4c< 2 -1] 
hence it is larger than N(N-2) for a > - 1. Noticing that both y* 2 /y 2 and || t || 2 K A / H^n+i Ho.k.a 
are bounded independently of N, we obtain the two inequalities of the lemma. 


The following result is proved in [BM 1, Lemma 2] in the particular but not simpler case 
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cx = 0. 

Corollary V.2 : For any real number a > -1 , the following inequalities hold for any 
polynomial v in P N _,(A) 

(V.26) o’’ N-' |( I -C 2 )v || 2 . A « ((1 -C 2 )v,( I -C 2 )v). 0 < c |( I -C 2 )v || 2 . A . 


Proof: Since v belongs to P N1 (A), we con write 

N + 1 

:2 ' -n 


w - Y" w n l* " 

v “ V 


By using twice formula (11.8), we know that 

( P«* 2 J n")" = (n-1)n(n.2ou 1 )(n*2ot«2) p R J* . 

Integrating twice by parts, we obtain 

(V.27) ||(1-e)v|| 0 % A ,L n N \! (v n ) 2 (n-1 )n(n + 2cx +1 )(n + 2a + 2) ||J“ ||q kA . 

On the other hand, noting that the Gauss quadrature formula is exact on P 2 n-i(A)' we have 
((1 -C 2 )v,( 1 -t. 2 )v) a0 =L n N ; 2 2 (v n ) 2 (n- 1 )n(n + 2a • 1 )(n+2a + 2) ||J* || 2 a A 
+ (v N ) 2 ((1-e)J N K " . (1-C 2 )Jn")«,o 

2\/~N- 1 ,-Ntl ,a 


(V.28) 


The two last terms have been computed in Lemmas V.2 and V.3 respectively. Consequently, we 
obtain the desired result by comparing (V.27) and (Y.28). 


We are now in a position to prove the 

Proposition V.3 : Let <x satisfy — 1 < oc < 1. For any integer £ , the form b* N satisfies the 
following properties of continuity 

(V 29) [ V w„ € P„, 3 (A) fl H 2 0 (A), V z„ € P„, 3 (A) D H 2 0 (A), 

W f. • 7 tl^ C 1 f ' ll'" N IIG n A 1 II ||j . A ) 2 ( ? l|z„ ||q . A 4 ||2 N ||| . A ) % ' 2 . 

and of ellipticity 


(V.30) 


Vw »« p K,3< A >nH 2 0 <A). 

b «,N< W N . W N ) > C ( f 4 N-’ !|w„|| 2 . A . P 2 I|W„||^. A , ||w N |||.. A ) 


Proof : Let w N and z N be two polynomials of P Nt3 (A) fl H a 0 (A). We study separately each of the 
three terms in b* N (w N , z N ). 

1) We have 

(w (lv) z ) - (w (lv) z ) 

VW N ’ ^n'u.,2 ~ ' W N > Z N ; «. ■ 
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so that Proposition 111.1 yields at once 
(V.31) (wjj ) , Z N ) a2 < c II W N ll2,ct,A II Z JI2,(X.A 


(V.32) (w^ * . w n ) k2 > c ||w N H 2 A 

2) We write the expansion 

w N = L n N : 0 3 W"J n ‘ and z^L^ 3 z"J n * ; 
since the quadrature formula is exact on P 2 n+ 3 ^^> we ^ ave 

- K • 2 N )«.2 ♦ K , Z N ) tt = - W N+3 Z N + 3 (J« + 3 , J N \ 3 ) K , 2 . 

Then, Corollary V. 1 yields on one hand that 

(Y.33) - ( W N ' W N^«,2 ^ “ ( W N ’ W N^« ^ C II W N 111 ,ot,A > 

and on the other hand that 

- < W N - Z N>*,2 = - K ■ Z N>« + C ^ N 2 ll J N + 2 llo 2 .a,A ' 

Noting that, due to (11.6), the coefficient of j£ +2 in the expansion of w N ' is equal to 
w N+3 (2N + 2oc + 5)(N+oc + 3)/(N + 2oc+3),we obtain 

- < W N * Z n)*, 2 < C < Kll2,«,A Mo,*.A + H W Ju,A ll Z N lll,oc.A > • 

Since H^(A) is the interpolation space of index 1 /2 between L 2 (A) and H 2 (A), this implies 

(V 34) - (w N . 2 n ) k2 ^ c (||w N || 2ct A ||z N || 0ocA + ||w N || 0 ot A ||w N || 2jKjA ||z n || 0 k A ||z n || 2jK A ) . 

3) Finally, writing w N = (1 -t; 2 ) 2 w N andz N = (1 -t 2 ) 2 z N . where w N and z N belong to P N _,(A), we 
note that, due to (II. 14), 

^ W N ’ Z N^et,2 = ^ 1 -C 2 ) 2 W N , < (( 1 ~^ 2 ) 2 W N . W n ) k+2 ,0 ^ 1 - ^ 2 ^ 2z N ' Z N ^«+2,0 

Applying Corollary V.2 with oc replaced by cx + 2 and noting that ||( 1 -t; 2 ) w N t| 0 a+2 A (resp. 
||( 1 -t; 2 )z N ll 0 ,et + 2.A :) coincide with ||( 1-t 2 ) 2 w N || 0etA (resp. ||( 1-t. 2 ) 2 z N || 0(t A ), we obtain 

(V.35) (w N . 2 N )«.2 ^ C il w N llo.ot.A ll Z N llo.*,A • 

On the other hand, Corollary V.2 implies 

(V.36) ( w n ’ w n^«, 2 ^ c N ||w n || 0 k A . 

■ Finally, the continuity property follows from (V.31), (V.34) and (V.35). The ellipticity 
property is a consequence of (V.32), (V.33) and (V.36). 

Theorem V. 1 : Let <x satisfy -1 cocci. For any f in C°(0) ,problem (V.18) has a 
unique solution u s in X $ . 
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Proof : Due to Proposition V.3, for any integer £ between -L and L, equation (V.20) has a unique 
solution in P N+3 (A) fl H 2 0 (A). The function u s □ £j._ l w^ exp(ifx) is the unique solution of 
problem (V. 1 8) in X 6 . 

We are going to prove an error estimate between u and u $ . If X is a Banach space, we 
introduce the projection operator n* from L 2 (0,X) onto $ L (0,X) which, with any function v with 
Fourier coefficients v*. t? € TL, associates the trigonometric series n'v = £<__ L v* exp(i£ x). We 
recall the following result [CQ, Thm 1.1 ], valid for any real numbers r and s, 0 < r ^ s : for any 
function v in H*(0,X), 

(V.37) IIv-n L v|| H y e X ) < c L II v H hj/o.x) 

A similar estimate [CQ, Thm 1.2] holds for the interpolation operator if s is > 1 /2 : 

(V.38) l|v-C L vlinyo,X) ^ c IMIhj(©,x) 


We conclude with the 

Theorem V.2 : Let a satisfy -1 <cx<1. If the solution u of problem (V.1) belongs to 
H* # (^) f° r a rea l number a ^ 2 , and if the data f is such that the function (1 -y 2 ) 3/2 f 
belongs to a space H* # (Q) for a real number q > 1 , the following error estimate between the 
solutions of problems (V.1) and (Y.18) is satisfied 

(V. 3 9) llu-u $ || 2 R# 0 ^c{(L 2 - ff 4 N 3 - ff )||u || ffa#0 4 N(L- f + N ,/2 " p ) ||( 1 -y 2 ) 3/2 f || # a # 0 } . 


Proof : We have 


(V.dO) Hu-u 6 || 2k#0 ^l|u-n'u || 2ci#0 . 

1) Using the definition (V.5) of the spaces H* # (Q), we have by (V.37) 


iiu'-a 


i2 

A 


)} 


1/2 


l|u-n L # u|| 2K#n < llu-n'u|| H j (0L g (A)) + \\u-n' u|| L 2 <e .H«(A)) 

< C L 2 ( ||u|| H j( 0 L g( A )) + H u HhJ-2(0,H^(A)) ) > 

whence 

(V.41) ||u-n*u|| 2itt(#(0 < c L 2 " ff IIulla.ot.^.rt ■ 

2) For any integer £ , -L ^ F < L, Proposition V.3, together with (V.9) and (V. 1 9), implies that, 
for any polynomial w^ in P N (A) fl H 2 0 (A), 
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I^N-WnIIo.oc.A + II'Jn-Wn Hf.ot.A < C N b a.N^N- W N ■ ^~ W n) 

< C N [bju'-wj, uj-wj) - (f\ Gj-wj). + d' f f , uj-wj) 2 ] 


(note that b* N (w^ * u^-w^) coincides with b*(w^ , u^-w^)), so that 

< ( 4 la'-ajig^ * iio'-o'in^ )•« < c n (f iiu'-w'ii 0 2 >a * iio'-w'ii 2 „ a )" 2 

, c- N sup — ,v ">« ,V " ) »' 2 ) 


*- c‘ N sup 

V N € ^N+3^^ 0 H^ 0 (A) 


'N "2,«.A 


To estimate the last term, we observe that 


^ > V n )« ^ > V n )«,2 ~ ^ “H ^ ’ V N^k + > V N^oi " (f* • V N ^«.2 

- «'<<’ . v H ). * (f'-HV. v N )„, 2 . 

so that, due to Lemma B.2, 

<?' • V* - ^ . v n )«. 2 < o {||( l -y 2 ) 3/2 «'-^')\\ 0 „ A ♦ N 1/2_p ||(i-y 2 ) 3/2 

+ N 1/2_p ||( 1-y 2 ) 3/2 A } II V N ll 2 ,ot,A . 

where we have chosen 1/2 < p‘ < p-1 /2. Choosing w^ = n^' 2 u ? and applying Theorem A.6, we 


derive 


(V.42) 


t -*112 \ 1 /2 


c N 3_<r (t? 4 ||II 2 


^-2,«,A * ) ,/2 + C-N(||(1 -y 2 ) 372 (?^-<f f )llo,.A 


4 N 172 - p ||( 1 -y 2 ) 372 ?'|| #ttA + N 1/2-p ||( 1 -y 2 ) 372 (f'-H # f*)|| p>A ) . 

The final estimate follows from (V.40), (V.41), (V.42) and (V.38). 


Remark V. 1 : By using the same techniques as in [BMM, § 2], one can check that, if the data f 
belongs to H p „(Q) for a real number p ^ 0,the solution u of problem (V. 1) is in H P+4 (Q), 


By analogy to Section IV, one could think of defining the following discrete problem 


Find u 6 in X 8 such that 


A 2 u s Cx k , t;f +,) = Kx k , -Ls$k^L, 1 ^ j < N 


(V.43) u 8 (x k ,±1) = 0 , -L ^ k < L , 

(0u $ /0n)(x k , + 1) = 0 , -L < k < L , 

(where the nodes are those of the Gauss-Lobatto formula). However, it turns out that this problem 
is not so accurate as (V. 1 7), Indeed, setting 

<P = Jjj/, - ((N*2a + 2)(Nioc)(N+<X4 1 )/(N-1 )N(N+2oc)) J“_; , 
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we observe that ((1 -t, 2 )<P ,tp) tt / Ho « +1 a’ s boun ded independently of N. On the other hand, 
for 1 < j < N, |«p(C*)| Is bounded by (c/N) |'(i;*)I. so that ((1 -C; 2 )«P ,<p)* >0 / Ho,« + i,a is 
^ c N' 2 . Noting that tp vanishes in + 1, we have found a polynomial ip = tp/( 1 -C 2 ) in P N1 (A) such 
that 

((1 -C 2 ) 3 *P .<p) tt0 < C N" 2 ((1 -C 2 ) 3 >p,<p) w 

Using this result with oc replaced by a +1, we obtain a polynomial w N in P n+ 3 (A) D H 2 0 (A) such 
that 

K > W n)«.1 < cN ' 2 H W nIIo 2 .«.A • 

Consequently, we have proved that the constant of el 1 ipticity of the bilinear form associated with 
problem (V.43) is s$ c N“ 2 . The convergence of the solution of this problem to the solution of 
problem (Y. 1) is not so good as the one we obtained in Theorem Y.2. 
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A ppend ix A ; Approximation erro r?. 

In this first appendix we are going to derive some results concerning the weighted Sobolev 
spaces we use in the paper and then we shall analyse the best possible fit in these spaces. The 
proofs we give generalize the analysis presented in [M] that focussed on the weight oc □ -1 /2 . 

A. 1. Results of interpolation between the weighted Sobolev spaces. 


We begin by stating some properties of the dual spaces H^ r (A) of H^ 0 (A). In the following, 
as is natural, we shall identify L*( A) with its dual space. As a consequence, the differentiation in 
the space of distributions 2>'{ A) is defined as follows: 

(A. 1) V f € 2>‘(A), V ip € £>(A), <df/dt, <p> n - < f, Q_ <Jt [d(<pp ot )/dC] > • 

Obviously the previous definition of the derivative coincides with the classical notion of 
differentiation for regular functions. 

Let us first introduce the space ^ r (A) of all derivatives of order r of functions of L*( A). 
The following theorem gives a characterization of the dual space H^ r (A) of H^ 0 (A), when L^(A) 
is identified with its dual space. It is well known in the case a = 0 and has been first established in 
the case oc = -1 /2 in [M, Thm 111.4]. 

Theorem A.1 : Let oc satisfy -1 < oc < 1. For any integer r, the spaces 3£‘ r (A) and 
H^ r (A) coincide. 


Proof : Let us first recall that, as a consequence of [BM2, Lemma III.2], the mapping: 
ip -> Q_ ot [d r (*p q^) /dC r ] is continuous from H^ 0 (A) into L*(A).Then, let us consider an elements 
of H^ 0 (A); it follows that for any f in L*( A), we have 
< f, Q_ ot [d r («pp ot )/dC r ] > <c||f|| 0(cA II*P ll r ,«, A • 

We derive from (A. 1) that #^ r (A) is contained in H^ r (A). 

Conversely, the range of H^ 0 (A) by the one-to-one mapping: <p -► p_ ct [d r (ipp K )/dt; r ] is a 
closed subspace of L*( A). Hence, any element L of the dual space H‘ r (A) of H^ 0 (A) defines on 
that closed subspace an element L such that 





- 42 - 


V<p€H; 0 (A), < L, ip > = < L, p_,jd r (tpp a )/dt; r ] > . 

Using now the Hahn-Banach theorem, we extend L to an element of the dual space of the entire 
space L*( A), which can be identified to a function f of L*( A) as previously mentioned. Using once 
more (A.l) allows to conclude that L and d r f/dC r coincide, which implies the imbedding 
H‘ r (A) c $£~ r (A). 


For any pair of Hilbert spaces X and Y such that X is contained in Y with a continuous 
imbedding and dense in Y, for any real number 0 in ]0,1 [, we denote by [X,Y] 0 the interpolation 
space of index 0 between X and Y, as defined in [LM, Chap. 1 ]. The following result of interpolation 
between the spaces H^ 0 (A) is proved in [BM2, Lemma III.3] (similar results can be found in [Gr] 
in a slightly different framework). 

Theorem A.2 : Let ex satisfy -1 < a < 1. For any real numbers p , q and s which satisfy 
0 ^ q < s ^ p and do not belong to IN + (I + cx)/2 , the following equality holds between the spaces 
of interpolation 

[H' 0 (A),H’ 0 (A)] (p ., |/(p . ql = H* 0 (A) . 

By duality, we derive the following result. 

Corollary A. 1 : Let cx satisfy -1 < a < 1. For any real numbers p , q and s which satisfy 
0 ^ q < s ^ p and do not belong to IN + (1 -*a)/2 , the following equality holds between the spaces 
of interpolation 

[h;«(a).h;'(a)] ( s .„ /Ib .„ =h;’(a) . 

As a consequence of [LM, Chop. 1, Prop. 2.1 ] we obtain 
Lemma A. 1 : Let ex satisfy -I < ex < 1. For any real number s which does not belong to 
IN i (Ua)/2 , the space [H* 0 (A), H‘*(A)] 1/2 is equal to L*(A). 

From this lemma it is simple to derive as in [LM, Chap. 1, Th. 6.2] the following general 
interpolation result. 

Theorem A.3 : Let cx satisfy -1 < oc < 1. For any positive real numbers p and q which do 
not belong to DM + (1+oc)/2 and for any 0 in ]0,1[,sef s = (1-0)p-0q. If neither s nor -s 
belongs to (N + (1 +oe)/2 , the following equality holds between the spaces of interpolation 


- 43 - 


[H« 0 (A), H“ q (A)] 0 = |h:,o(A) if $> 0 , 

H*(A) if s <0 . 

Now, we want to extend this result to the entire weighted Sobolev spaces. 

Theorem A.4 : Let oc satisfy -1 < <x < 1. For any real numbers p and q which satisfy 
0 q < p and for any 9 in ]0,1[ , the following equality holds between the spaces of 
interpolation 

[Hj(A), H;<A)] 9 = H^ 1_0)p+0q (A) . 

Proof : a) We begin with the case where q Is equal to 0 and p and 9p are integers. We first note 
that the operator d p /dt; p is linear continuous from H P (A) into L*(A), besides, from Theorem A. 1, 
it is also linear continuous from L*(A) into H^ P (A). The principal theorem of interpolation 
[LM.Chap. 1 , Th. 5.1] states that d p /dt p is linear continuous from [H P (A), L*(A)] 0 into 
[L*(A), H‘ p (A)] e . From the previous theorem, this space coincides with H‘ 0p (A) so that we 
obtain the imbedding 

(A.2) [H P (A),L2(A)] 0 c{f€L^(A) ; d p f/dt; p €H; 0p (A)} . 

Using now the characterization of the space H^ 0p (A) given in Theorem A. 1, we derive that, for any 
f in L*(A) such that d p f/dC p belongs to H^ 0p (A), there exists an element g in L*(A) such that 
d p f/dt p =d 0p g/dt; 0p . This means that 

d 0p ( (d ( 1 _0)p f/dt; ( 1 _e)p )-g)/dt; ep = 0 , 

hence that (d (,_0>p f/dt; n “ 0>p )-g is a polynomial of degree < 0p- 1, thus an element of L*(A). We 
conclude that 1 - e, Pf/dt; < 1 _0)p is an element of L*(A) ; hence, from (A.2), we deduce 
[H P (A), L^(A)] 0 C { f € L*(A) ; d (1 - 0)p f/dC (1 " 0)p € L*(A) } = H' 1_0)p (A) . 

Let us recall now that the reverse imbedding was deduced in [BM2, Lemma 111.3] from the results 
of Theorem A.2. This proves the theorem in this special case. 

b) The general'case is derived by using the reiteration theorem [LM, Chap. 1, Th. 6.2] in 
three steps. First, for any integer p and any real number 0 in ]0,1[, denoting by m the integer 
such that m- 1 < (1 -0)p < m < p, we have 

H^A) = [H;<A),Hr'CA)] m . M . e)| , 

= [ [H:<A>, L 2 (A)],. m/p , [h;(A), L 2 (A)],_ {m _,, /p = [Hj(A),L 2 (A)] e . 

Next, for any positive real number p and any real number 9 in ]0,1 [, denoting by s the integer 
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such that s- 1 < p < s, we write 

[h:<A),L„ 2 (A )] 9 = [ [H;(A),L 2 (A)],. p/> ,[h:(A),L 2 (A)], ], 

= [H:(A).L 2 (A)],.„. e)p/1 =Hi , - , »'(A) 

Finally, for any positive real numbers p and q, p ^ q, we derive 

[h;<a), h;<a)], = [ !h;(a),l 2 <a)j 0 , [hS(a),l 2 (a)],.,, /p ], 

= tH'(A),L. 2 (A)], ( ,_ p/p) = H:'- M - 9 «(A) . 

It follows from [BM2, Lemma 11.5] that, for any nonnegative real number r which does not 
belong to IN + (1 foO/2, H^ 0 (A) is the space of all functions in H^(A) which vanish in ± 1 
together with their derivatives of order < r - (1 -toO/2. Thus, the following corollary is a simple 
application of [LM, Th. 13.3]. 

Corollaru A.2 : Let cx satisfy -1 <<x< 1. For any real numbers p , q and r which satisfy 
0 < r < q ^ p and such that r does not belong to IN + (1 +oc)/2 and for any 0 in ]0,1 [ , the 
following equality holds between the spaces of interpolation 

[h£(a) n h; 0 (a), h;(a) n H; 0 (A)] e = H^- e>p+ 6 q (A) n h; 0 ca) . 

We are now in a position to study the best approximation errors in several norms. 


A.2. Analusis of the best approximation in weighted norms. 


We begin by analysing the best approximation of elements of H^ 0 (A) by polynomials of 
P n (A) fl H^ 0 (A), when r is an integer. This is the most simple situation and can be done by 
induction. 

Theorem A.5 : Let ex satisfy -1 < <x < 1. For any nonnegative integer r , the orthogonal 
projection operator n£ ^ from H^ 0 (A) onto P N (A) fl H^ 0 (A) is such that , for any real 
number a > r and for any ip in H"(A) fl H^ 0 (A) , the following inequality holds 
(A3> ll» —<>IL A « cNr "’IM,,«. A • 


Proof: For r = 0, the result is well known and has been proven in [CQ, Thms 2.1 and 2.3] for the 
cases oc = 0 and a = - 1 /2 and in [BM2, Thm IV. 1 ] for the general case. The following will be an 
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induction process as we suppose that (A.3) holds for any integer N, with r replaced by r-1 and for 
any real number o > r-1. Then, for any function <p in H*(A) fl H^ 0 (A), we note that <p' belongs 
to H^'Va) f) H^"q(A). Let us define ip N as 

<P N (0 = f_\ [nS:^. 1 (<p')(0-m(nJ:;- r ^/<p1)(1-L 2 ) r ‘ , /m((1-t; 2 ) r ' , )]dL , 

o r 1 

where, for any ip in L K (A), m(ip) stands for the integral J_, ip(L) dt;. It is simple to note that <p N 
is a polynomial of degree ^ N. Besides it vanishes at ± 1 , as well as its first (r- I) derivatives. 
Using now the Poincare-Friedrichs inequality, which is valid in the weighted Sobolev spaces as is 
proven in [BM2, Corollary III. 1 ]: 

we deduce that 


I1Hr,«,A < c (||(p - n N ^ t r _,(<p)^ A + | m(n N ^ r _,(ip)) |) 

<oik-nSj; r , . 1 (f) \\ r . UA ■ 

Using now the induction hypothesis, we deduce that 

II* - *H l,.A < 0 N''-"-'-" Ilf ||„_, , A < 0 N'- U ll,„. A . 
which complete the induction argument since ||(p - njj r<p|| raA ’ s defined as the minimum of 
II <P - <P N ll r>K A over all the «p N in P N (A) fl H* 0 (A). 


Remark A. 1 : It is Important to note that the operator does not have optimal approximation 
properties in higher norms than II ll ret A • For instance it is proven in [CQ, §2] and in [BM2, § IV ] 
that there exists a function tp in H^fA) such that the L 2 -projection operator satisfies 
ll^-~<0%ll,. tt .A>oN 1/2 |kll 1 ., A . 

which legitimates the introduction of orthogonal projection operators in H^ 0 (A) for any integer 
r ^ 0. One important question still remains : what are the approximation properties of the 
operator ^ in lower order norms? The following theorem states that it remains optimal. 


Up to now, we have not explicitly used the exact formulation of the norm on H^ 0 (A) to 
define the operator . In order to be able to perform a duality argument, we shall write it. As a 
special case of Theorem A.2, H^ 0 (A) can be seen as the space of interpolation of index 1/2 
between L 2 (A) and H 2r 0 (A). Consequently, there exists an unbounded self-adjoint linear 
operator ©, which is positive definite and such that the domain D(0) of 0 in L 2 (A) is H^ 0 (A) 
and the domain D(0 2 ) of 0 2 in L 2 (A) is H 2r 0 (A). This allows for defining a scalar product ((.,.)) 
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cver H^ 0 (A) as follows : 

(A.4) ((u,v)) = (u,v) a + (0u,Ov) R . 

The corresponding orthogonal projection operator will be denoted by . The interest of this 
new scalar product is related to the regularity property of the following operator : let g be in 
L^( A), Tg is the solution of the problem 

Vv€H «,o (A) > (( t 9. v )) = (9 i v ) 0 ,«,a • 

Then it is simple to check that T is linear and continuous from L*(A) into H^ r 0 (A). This is the 
main ingredient of the duality argument to derive an estimate of tp — tc« p in the L^( A)-norm. The 
reader is refered for instance to the abstract duality theorem of [M, Thm A. 1 ] for more details. 
The scale of all estimates of f-njjf in the H*(A)-norm, 0 < s < r, is then derived by 
interpolation and we obtain the following theorem. 


Theorem A.6 : Let a satisfy -1 < oc < 1. For any nonnegative integer r , the orthogonal 
projection operator n^ r r from H^ 0 (A) onto P N (A) f] H^ 0 (A) for the scalar product defined 
in (A.4) is such that , for any real numbers a and s , 0 s < r < a , and for any f in 
H"( A) H H^ 0 (A) , the following inequality holds 
(A.5) ||f - njj ptp \\ 9a A < c N ,_<T ||f ||„ k A • 


Remark A.2 : Theorem A.6 is valid for any positive real value of r as can be derived from the 
same arguments as in [M]. We shall not present this extension which is not used in this paper. 


The next step is the analysis of the best approximation of elements of H^( A) D H£ 0 ( A) by 
polynomials of P N (A) Pi H£ 0 (A), when r and p are integers, p < r. Let f be an element of 
H"(A) ft HJ 0 (A) and f 0 the polynomial of P 2r _,(A) such I 1 ' fo be t° n 9 s to H^ 0 (A). It is 
simple to note that tt^ r(<P-<P 0 ) + <P 0 is an optimal approximation of f, more precisely, 


, II«p -«; r r (f-<c 0 ) + t P 0 ]ll s .«,A < cN 

This leads to 


3-<r 


IV a, A 


<c’ N 


S-<J 


<r,ot,A 


Theorem A.7 : Let a satisfy -1 < a < 1. For any nonnegative integers r and p r , there 
exists a projection operator from H^(A) fl H£ 0 (A) onto P N (A) 0 H£ 0 (A) such that , 

for any real number o > r and for any f in H®(A) f) H£ 0 (A) , the following inequality holds 
(A.6) ||f - njjf || raA < c N r ‘ 9 ||f ||, KA 
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Remark A.3 : Here again, the result can be extended to any real value of r as in [N]. 

Finally, we can state the 

Corollaru A.3 : Let oc satisfy - 1 < a < 1 , For any nonnegative integers r and p , for any 
real number a^r and for any <p such that the function (1-t; 2 ) p <p belongs to H"(A) , there 
exists a polynomial ip N of P N (A) such that 
(A.7) || (l-t 2 ) p (ip-ip N )ll r .«.A < 0 N r -' II (l-C S ) p II, , A . 

Proof: The result is completely obvious and results from Theorem A.7 in the case where p is < cr 
since <p N can be chosen equal to [n^;^ j(( 1 -C 2 ) p <p)]/( 1 -C 2 ) p . where a denotes the integral part 
of o. Otherwise, we define (1 -t; 2 ) p <P N as the orthogonal projection of tp onto P N+2p ^^ fi H£ 0 ( A) 
for the scalar product associated with the norm : <p -► || (1 -t; 2 ) p <p ll r w A and we obtain the result by 
interpolation between the two estimates 

II (1 -t 2 ) p (f->P N )|l r .«. A «II Cl-C 2 ) D >p|| r ... A 

II ('-C 2 )"(»-* N )ll r .,. A < ll<l-t 2 )%-[n N X<C-t J >N>]lU A 

< c N r -" ||( 1 -C 2 )MI P ,. A ■ 


and 




Appendix B ; Interpolation errors. 


This appendix deals with the approximation properties of the quadrature formula analysed 
in the first section. Let us introduce the related operator of interpolation at the internal 
points, defined for any function <p continuous on A as follows ; belongs to P N _,(A) and 

satisfies 

(B.l) . t<j<N . 

The following estimate is well-known in the case a = 0 [CQ, Thm 3.2], it can easily be proved for 
any ex > -1 in a similar way : for any real number a > sup {1 /2, (1 tcO/2} and for any function 
<p in H"(A) 

(B . 2) || <p _, <P || 0 K A < c N 1 /2_(r || <p || ff K A . 


The previous inequality implies as a special case that, for any integer m > 1 and for any 
smooth function <p, the quantity ||( 1 -C: 2 ) m/2 (tp —) Ho«a (ends *° 0- following lemma 
gives a precise form of this result. 

Lemma B.l : Let a satisfy -1 <a < 1. For any integer m 1 , define k as being equal to 
m/2 or (m + 1)/2. For any real number a > sup {1 /2, (1 +oO/2) and for any <p such that the 
function (I -L 2 ) k <p belongs to H"(A) , the following inequality holds 
(B.3) ll(l-C 2 ) k (»-t;.74i)|| 0 , A <CN' /2 -’||(l-!; 2 )%ll,,,, A . 


Proof: Let <p N be any element of P N .,(A). We use the triangular inequality 

<B.4) IK I -C 2 ) k s IK l-c 2 ) k <«-«„) ll„,«, A ■ II (i-t; 2 ) 1 ll 0 ,. A . 

Let us recall that the quadrature formula (.,.) Km is exact on * n fh ecase 2k = m, we 

deduce that 

ii< i-c 2 ) k <♦„- o) ii 0 2 ... a = IJ 1 ., i-(cp m ) 2 ) 2t . 

In ttys case 2k = m + I, we have 


IK 1-C ) f‘P n - ^N_ 1 llo,«,A -IKl-C ) (<Pn ^N -1 ( P'" 0 .«+m,A • 

Using exactly the same techniques as for [BM1, Lemma 1 ], we prove that, for any <x > -1 and for 


.2x1/2 
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z:,"„ (i-<t,*) 2 ) ♦„*«;;) oj 


.oc,0 


= L" n :\ cx 2 n(n + 2 « + 1) \UX«A + < £j\i <Mt*) 2 ) J^CtJ 1 ) 
□ Ll'! cx 2 n(n + 2a + 1) ||J“|| 0 2 aA + a 2 N(2N + 2cx + 1) 

> ZL « n 2 n(n + 2cx+1) ||J*IIq K A = ll(1-t 2 ) ,/2 4- N ll 0 2 . tt . A • 


This estimate, with cx replaced by cx 4 m and >j/ N = <p N - ■ 9' ves 


-e> k (♦«- < l-,, ( 4 - 4 > N ) 2 (cr m ) ,r° 


^ V N /« . \2/>-«+m\/, / y «+m\2\2k -«,m 

< L j= i (tp-4> N ) (Cj KMtj ) ) Pj 


Consequently, in both cases, we obtain 

-2\k / . . oc+m ii2 


n< i -c 2 ) 1 <f H - O) ii 2 A < in i -c 2 ) 1 **' 2 <♦-♦„•> ii, 2 -, Al z-,, < i -(t;* m > 2 r p; 

If ex is < 0, it is an easy matter to prove that the sum on the right-hand side is bounded : indeed, 
all the (1 -(C**™) 2 )'* are < 1 and all the corresponding weights are positive; since the constant 
polynomial is exactly integrated by the formula, we derive 

L N / , «+m\2s-« „«,m - v ^ ~«,m o 

j=1 (1 “(Cj ) ) Pj ^ L js i Pj < 2. 

If a satisfies 0 ^ cx < 1, we recall that, due to (11.14), we can write 

z“._, < i -(c,- m ) 2 )-* p,' m =Lj"., (i -«r) 2 >’" p,-'' m4 ' 

and the same argument as before can be used since now 1 -cx is ^ 0. In conclusion, we can state that 
there exists a constant C such that 

ll( 1 -C 2 ) k <♦„- II?. „.a < 0 IK I -t 2 > k,,/2 <♦-♦„> »L~<AI ' 

The Gagliardo-Nirenberg inequality now yields to the bound 

IK I -t 2 ) k <»H- O llo.«A « C IK I-i: 2 ) 1 *" 22 (♦-♦„) lj! W IK I-i; 2 ) k * </2 <♦-*„) llnl(A) ■ 

Besides, since the multiplication by p_ w/2 is an isomorphism from L 2 (A) onto L 2 (A) and from 
Hq(A) onto 0 (A) [BM2, Lemma 111.2], we derive that 

(B.5) ||(i-t 2 ) k (<P N - ^_7<p) ||J iKiA < C||(i- t 2 ) k (< p-«p n )II 0 ,«,a IK 1 -t 2 ) k («P-<P N ) II 1a , a • 

The inequality (B.3) is now a simple consequence of (B.4), (B.5) and Corollary A. 1. 


- oc.-hti\ 2\-(X * <x,m 


We can state the properties of the quadrature formula. 

Lemma B.2 : Let cx satisfy - 1 < a < 1. For any integer m > 1 and for any real number 
a > sup (1 /2, (1 >oc)/2} , for any f such that the function (1-C 2 ) <m+1>/2 f belongs to H"(A) , 
the following inequality holds 


(B.6) sup 

V N € ^N+ 2 m- 1 ^^ ^ H Ri 0 (A) 


^■ V N^» ~ ^ ’ V N ) Km < c | s | 1/2 - <r 


llv N ll 


_^2ym +1 )/2 


fll 


vaA ‘ 


N Hm.oc.A 
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Proof : Since the quadrature formula is exact on P 2N+2m _ we have for any v N in 
P N + 2 m-i (A) H H k " 0 (A) and any f N in P N _,(A) 

^> V N^oc “ ^’ V N^«.m = ( f-f N ,V N^« “ ^N-1 ,V N^oc 

Next, we note that v N can be written (1 -C 2 ) m v N , where v N belongs to P N-1 (A). Moreover, using 
once more the fact that the multiplication by q^ /2 is an isomorphism from H^ 0 (A) onto Hq(A) for 
any k ^ 0 and that the multiplication by (1 -C 2 )" m is continuous from H™(A) into L 2 (A) [LM, Th. 
11.2 and 1 1.3], we have 

II^N llo,«,A = II^N Pot/2 llo.A ^ C IK 1 Pot/2 Ilm.A = C ll V N Poc/2 Ilm.A ^ C ll V N Ilm.oc.A 

Thus, we derive 

(f- v N ) a - (i’.V N ) am = (f-f N ,V N ) a+m - AN^ot+m 

^ C [IK 1 llo,cc,A + 11^ ^ ^N- 1 IIo,«,a1 II^N llo.ot.A 

<c[||(1-C ) (f-f|p llo,«,A + IK i ) (Hi- 1 llo.ot.A^ ll V N IL,ot,A 

Since (1 -C 2 ) m is upbounded by (1 -t; 2 ) <m+1)/2 , we obtain 

(f,V N )*-(f,V N ) K ,m 

<o[||(i-t 2 ) lm>,,/2 (f-f„)|| 0 .. A *|(i-c 2 )'"*' ,,2 (tj : ” f -g II 0 .«. a ]IIv„IU a . 
and the result follows from Lemma B. 1 together with Corollary A. 1. 

Remark B. 1 : We actually proved the stronger result 

sup <cN ,/2 -’iKi-c 2 ) lm, ” /2 fiu. A . 

V N € P N+2m-1 ( A ^ H H K<0 (A) V N ISf.ot,A 

where £ is equal to (m ♦ 1 )/2 if m is odd and is > (rm 1 )/2 if m is even. However, we do not need 
this improvement. 


Remark B.2 : in the previous analysis, we only need the Lagrange interpolation operator. With 
the generalized Gauss type formula, it seems natural to associate the Hermite interpolation 
operator 3^ +2m _ t : for any function <p in C m_, (A) I belongs to P N+2m _i( A ) ancl satires 

(B.6) (3jJ + 2 m _i<p)(C* +m ) = <p(L| t+m ) , 1<Um , 

(d k /dt; k )(3“ +2m _ 1 .p)(± 1 ) = (d%/dt; k )(±l) , CUk^m-1 . 

However, the approximation properties of this operator seems poor when m is larger than 1. 
Indeed, in the case m = 1, the following estimate [BM2, Lemma V.9] holds for any function ip in 
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H*(A), o > sup {1/2, (1 4 oc)/2}, 

(B.7) ll'P-(3I.,'f )II 0 ,».A « 0 N ' /2 ’’ M„.„,a 

But in the cose m = 2, due to Corollary V.2, the only interpolation result we are able to prove by 
using the same techniques is that, for any function ip in H"(A), a > sup {3/2, (3+a)/2}, 

(B.8) llf-(3;.,H)ll 0 ,,A«<=N ,/2 -"|kll,,A 
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Appendix C ; Nodes and weights. 

We end this paper by giving numerical examples of the nodes and weights involved in the 
generalized quadrature formula (II. 10). Let us briefly recall how they are computed (we refer to 
[GK][KE] for more general techniques). First, setting for a while J“* = J* / ||J“|| 0 kA , n <E (N, we 
have by (11.4) and (11.5) 





o o ... o 




0 




e, o e 2 ... o 


J- 


0 

(C.1) K 

... 

= 

o e 2 0 ... 0 


... 

+ &N 

0 


^.r 


° . 6n-iO. 


,4Li* 


l«* 

J N J 


with 


1 / n(n+2a) 

(C.2) &. = — - and & = / - - - , n > 2 . 

7 200 3 7 4 (n+oc) 2 - 1 

Hence, the nodes L* , I < j ^ N, of the Gauss formula are simply the eigenvalues of the previous 
symmetric tridiagonal matrix; due to formula (11.1 2), the weights p"'° , 1 < j < N, are equal to 
H^o Ho* a Xji 2 ’ w * iere Xj, is the first component of one of the corresponding eigenvectors with 
euclidean norm equal to 1. Therefore, nodes and weights can be computed for instance by a standard 
QR-algorithm (see [WR]). 

The internal nodes L“ +m and weights p“ ,m , 1 < j < N, of the generalized quadrature formula 
(11.10) are obtained by using exactly the same algorithm with a replaced by a + m, then applying 
formula (11.14). As far as the boundary weights are concerned, the first pair p£ - ] is calculated 
from formula (11.20). Then, for m ^ 2, the other pairs p*'™ , 1 < k < m-1, are obtained by 
induction on m, from the triangular linear system (11.21)(11.22); finally, the pair p^'™ is 
computed from formula (11.29) or (11.31). 

i 

Figures 1,2 and 3 represent the nonnegative zeros of the polynomial Jjj respectively for N 
= 6, 9 and 12, as a function of a, - 1 < a < 3 (recall that we are mainly interested with values of 
a between -1 and 1, and values of m equal to 0, 1 or 2, so that the range of a + m is ]-l ,3[). It can 
be observed that, when a grows from - 1 to 3, these zeros decrease slowly (but not linearly). 
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Figures 4, 5 and 6 represent all the zeros of the same polynomials for several values of a, 

namely for a = -3/4, — 1 /2, - 1 /4.11 /4, 3. As is well-known, for fixed values of oc and N, 

these zeros are not evenly distributed, but they cluster in the neighbourhood of ± 1. 

Finally, Tables I to X can be used to compute , $(0 P K (0 dL, for oc = -1/2, -1/4, 0, 

1 /4 and 1 /2, by a quadrature formula involving respectively 1 2 or 40 internal nodes; each table 
gives the nodes and weights of formula (II. 10) for m = 0, 1 and 2. Note that, as foreseen, the nodes 
of the formula for a = - 1/2 and m = 1 or 2 coincide respectively with those of the formula for oc 
□ 1/2 and m = 0 or 1. As it is well-known, in the Chebyshev case oc = -1 /2, for m = 0, the 
internal weights are all equal to the same quotient of it; the same is true in the case m = 1 and 
moreover the boundary weights are the half of the internal weights; these properties are no longer 
valid when m is equal to 2. 

Example C. 1 : Let us approximate the quantity 3 = (n/4) cos (nt;/2) dC = 1 by using the 
quadrature formula (II. 10) with N = 6 internal nodes, successively for m = 0, 1 and 2. We obtain 
m □ 0 : 3 « 0.9999999997386354 
m a 1 : 3 a* 0.9999999999989776 
rn = 2 : 3 « 0.9999999999999963 . 
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The six nonnegative zeros of J? 2 as a function of ex. 
















labia l 

The nodes and weights of three quadrature formulas with N = 12 internal nodes for a = -1 /2. 
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Table IV 

The nodes and weights of three quadrature formulas with N = 12 internal nodes for a = 1 /4. 
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Table V 

The nodes and weights of three quadrature formulas with N = 12 internal nodes for a = 1 /2. 
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Generalized Gauss formula (m = 2) 
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Table X 

The nodes and weights of three quadrature formulas with N = 40 internal nodes for a. = 1 /2. 
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